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ABSTRACT 

The absence of radio pulsars with long periods has lead to the popular notion of a high P "death line." 
In the standard picture, beyond this boundary, pulsars with low spin rates cannot accelerate particles 
above the stellar surface to high enough energies to initiate pair cascades, and the pair creation needed 
for radio emission is strongly suppressed. In this paper we explore the possibility of another pulsar "death 
line" in the context of polar cap models, corresponding to high magnetic fields B in the upper portion of 
the period-period derivative diagram, a domain where few radio pulsars are observed. The origin of this 
high B boundary, which may occur when B becomes comparable to or exceeds B CI = 4.4 x 10 13 Gauss, 
is also due to the suppression of magnetic pair creation, but primarily because of ineffective competition 
with magnetic photon splitting. Threshold pair creation also plays a prominent role in the suppression 
of cascades. We present Monte Carlo calculations of the pair yields in photon splitting/pair cascades 
which show that, in the absence of scattering effects, pair production is effectively suppressed, but only 
if all three modes of photon splitting allowed by QED are operating in high fields. This paper describes 
the probable shape and position of the new "death line," above which pulsars are expected to be radio 
quiet, but perhaps still X-ray and gamma-ray bright. The hypothesized existence of radio-quiet sources 
finds dramatic support in the recent discovery of ultra-strong fields in Soft Gamma-ray Repeaters and 
Anomalous X-ray Pulsars. Guidelines for moderate to high B pulsar searches at radio wavelengths and 
also in the soft and hard gamma-ray bands are presented. 

Subject headings: gamma rays: theory — radiation mechanisms — magnetic fields — stars: neutron — 
pulsars: general 



1. INTRODUCTION 



O Most current theories for the generation of coherent radio 
^eSnission in pulsar magnetospheres (see summations in Michel 
5©91; Melrose 1993; Lyutikov et al. 1999) require formation of 
.aii electron-positron pair plasma. Such a pair plasma has been 
Jghown to develop via electromagnetic cascades along the open 
agnetic field lines as a result of particle acceleration to TeV 
ergies, followed by curvature, synchrotron radiation and one- 
Jjroton pair production in the strong magnetic field, 7 — * e + e~ , 
near the neutron star surface (e.g. Sturrock 1971; Ruderman & 
Sutherland 1975, Arons & Scharlemann 1979). Numerical sim- 
ulations of pair cascades above a pulsar polar cap show that 
10 3 — 10 4 pairs are produced for each accelerated primary elec- 
tron (Daugherty & Harding 1982) in a typical young pulsar. 
But if any of the necessary conditions for development of a pair 
cascade, i. e. production of 7-ray photons (requiring particle ac- 
celeration to high energy) or production of sufficient pairs, is ab- 
sent the pulsar should not, according to theory, emit detectable 
radio emission. The existence of an observed "death line" in 
the P — P distribution of known radio pulsars (see Figure |5| 
below), to the right of which no pulsars were detected, until the 
recent discovery (Young, Manchester & Johnston 1999) of the 
8.5 second pulsar PSR J2144-3933, provided some circumstan- 
tial evidence that pairs are required for radio emission. The 
slope of this putative boundary fits a line of constant voltage 
across the open field lines, V oc p- 3 / 2 p 1 / 2 ( suggesting that 
radio emission ceases when the particle acceleration is not high 
enough to sustain pair production. The radio emission from 
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some long period pulsars like PSR J2144-3933 can be explained 
through production of pairs by photons from inverse- Compton 
scattering (Arons 2000, Zhang, Harding & Muslimov 2000), a 
process which does not require as high a voltage as does curva- 
ture radiation. 

In this paper, we study another area of pulsar phase space, 
that of very high magnetic fields, where the character of pair 
cascades is significantly altered by the process of photon split- 
ting and other effects such as ground state pair creation and 
positronium formation. Magnetic photon splitting, 7 — > 77 , 
a QED process in which a single photon splits into two lower- 
energy photons (Adler 1971, Baring & Harding 1997), operates 
efficiently and competes effectively with pulsar pair production 
only in magnetic fields above ~ 10 13 Gauss (Harding, Baring 
& Gonthier 1997, hereafter HBG97). This region of high mag- 
netic field strength lies in the upper-right part of the P — P 
diagram. We will explore whether photon splitting may be the 
explanation for why there are no classical radio pulsars detected 
with derived magnetic fields above ~ 10 14 Gauss. Our initial 
results (Baring & Harding 1998) have suggested that photon 
splitting could indeed explain the absence of very high-field ra- 
dio pulsars. This paper considers in more detail the underlying 
assumptions made in that calculation, explores the physics of 
high-field pair cascades and makes predictions for high-energy 
searches for high-field pulsars. Although pulsars in this region 
of phase space may be radio-quiet, those that are to the left of 
the death line are still prolifically accelerating particles to high 
energies. Their pulsed emission may therefore be detectable in 
the X-ray and gamma-ray wavebands, making this relatively- 
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unexplored region of pulsar phase space extremely interesting. 

There has recently been growing evidence for another class 
of pulsars with ultra-strong magnetic fields. Observations at 
X-ray energies have yielded detections of both long periods and 
high period derivatives in two types of sources, anomalous X-ray 
pulsars (AXPs) and soft 7-ray repeaters (SGRs), which suggest 
dipole spin-down fields in the range 10 14 - 10 15 Gauss. The 
AXPs are a group of seven or eight pulsating X-ray sources with 
periods in the range 6-12 seconds, and are continuously spinning 
down (Vasisht & Gotthelf 1997). The SGRs are a type of 7-ray 
transient source that undergoes repeated bursts; four (possibly 
now five: see Cline et al. 2000 for SGR 1801-23) are currently 
known to exist. A detailed discussion of these sources and re- 
cent discoveries pertaining to them is presented in Section [l.3j . 
While conventional radio pulsars tap their rotational energy to 
power their emission, if such a new class of ultra-magnetized 
neutron stars or "magnetars" does exist, their long periods in- 
dicate that their emission cannot be rotationally-powered: they 
may instead be driven by their magnetic energy. With the pos- 
sible exception of the recent report (Shitov 1999; Shitov, Pu- 
gachev & Kutuzov 2000) of a low frequency pulsed detection of 
a counterpart (PSR J1907+0919 at 111 MHz) to the soft gamma 
repeater SGR 1900+14, which has not been confirmed by pulsa- 
tion searches in Aricebo data at higher frequencies (Lorimer & 
Xilouris 2000) none of these sources has detectable pulsed radio 
emission. 

In several previous papers, we have computed the photon 
splitting attenuation lengths of photons emitted parallel, or at 
small angles, to the magnetic field near the surface of strongly 
magnetized neutron stars and compared them to the attenua- 
tion lengths for one-photon pair production (Harding, Baring 
& Gonthier 1996; 1997). Both processes depend on the photon 
energy w and the angle #kB of photon propagation to the local 
field. However, since photon splitting has no threshold, it may 
occur before the photon reaches the threshold for pair produc- 
tion at to sin 6*kB = 2m e c 2 . We found that the photon splitting 
attenuation lengths are lower than those for pair production in 
fields higher than ~ 3 x 10 13 Gauss, for emission in the open field 
region near the neutron star surface. One must therefore incor- 
porate photon splitting in cascade models of highly-magnetized 
pulsars. This has been done in the case of PSR1509-58 (Hard- 
ing, Baring & Gonthier 1997), which has a field of 3 x 10 13 Gauss 
as estimated using Equation (||). We considered both the case 
where all linear polarization modes of photon splitting allowed 
by CP invariance operate, and the case where only the _L— >|||| 
mode (our polarization convention is defined near the begin- 
ning of Section |J) allowed by selection rules (Adler 1971) in the 
weakly dispersive limit operates. The latter case, we argued, 
was most appropriate for pulsars. In both cases, we found that 
photon splitting attenuation is a plausible explanation for the 
very low spectral cutoff between 10 and 30 MeV in PSR1509- 
58 (Kuiper et al. 2000). Since the rates of these attenuation 
mechanisms are very sensitive to the polarization states of the 
incoming photons, it is evident that the polarization properties 
of the radiated photons are very important in modeling both 
the output spectrum and pair suppression in high-field pulsars. 

In ultra-strong magnetic fields, the vacuum dispersion in- 
creases to the point where the selection rules derived for weak 
linear dispersion (appropriate for B <J B cr ) may no longer be 
valid, and higher order non-linear contributions to such disper- 
sion may be significant. We presently do not know which modes 
of photon splitting operate in supercritical fields and thus do 
not know which photon polarization modes are allowed to split. 



We will therefore explore several possibilities for photon split- 
ting mode behavior in our study of high- magnetic field pair cas- 
cades. In Section 2.1, we discuss the physics of pair suppression 
in high magnetic fields, both by photon splitting and by pair 
production in low- lying Landau states, as well as the possible 
alternative mechanism of pair suppression by bound-state pair 
creation, discussed by Usov & Melrose (1996). In Section ||, we 
present results of detailed simulations of splitting/pair cascades 
which explore the conditions for suppression of pair creation 
in a pulsar magnetosphere. We compute the boundary in the 
pulsar P — P distribution where the escape energies for pho- 
ton splitting and magnetic pair production for photons of _L 
polarization are equal. This condition then defines a putative 
radio-quiescence boundary, posited in Baring & Harding (1998), 
above which pulsars may not produce the dense pair plasmas 
required for radio emission. For reasonable assumptions about 
the location of particle acceleration and the angles of the emit- 
ted photons in high-field pulsars, the computed radio-quiescence 
boundary lies at P above those in the known radio pulsar pop- 
ulation, potentially explaining the absence of radio pulsars with 
fields above ~ 10 14 Gauss. Our cascade simulations indicate 
that such a boundary is viable only if photons of both polariza- 
tions can split; otherwise pair creation is generally postponed a 
generation rather than suppressed. In the latter case, significant 
reductions in pair production can ensue only if the maximum 
energy of primary photons is below or not much above the pair 
creation escape energy. If pair suppression is rife at these high 
P , it suggests the possible existence of a new class of radio-quiet, 
ultra-magnetized pulsars that may be related to the emerging 
class of magnetars. 

2. PHYSICAL EFFECTS INFLUENCING PAIR CREATION 

The most important competitor to pair creation 7 — > e + e~ 
at high magnetic field strengths as a mechanism for attenuating 
photons in pulsar magnetospheres is magnetic photon splitting 
7 ~~* 77 ■ Hence this process motivated the suggestion (Bar- 
ing & Harding 1998) that splitting can potentially suppress pair 
creation and inhibit radio emission in highly-magnetized pul- 
sars; accordingly it forms the centerpiece of the discussion of 
this Section. Yet, the creation of pairs in low Landau levels or 
the ground state can also suppress sub sequent pair creation, so 
this effect is addressed in Section 2.2. We also include a brief 



discussion of the role of positronium formation, since the res- 
idence of pairs in neutral bound states can inhibit a range of 
coherent mechanisms for producing radio emission. The section 
concludes with a discussion of how the combined effects of pho- 
ton splitting and threshold pair creation change the nature of 
pair cascades as the magnetic field increases. 

At this point, it is appropriate to identify conventions adopted 
throughout this paper. The photon linear polarizations are such 
that || refers to the state with the photon's electric field vector 
parallel to the plane containing the magnetic field and the pho- 
ton's momentum vector, while _L denotes the photon's electric 
field vector being normal to this plane. Furthermore, all photon 
and electron energies will be written in dimensionless form, be- 
ing scaled by m e c 2 . Since general relativistic effects will play an 
important role in our considerations, we adopt the convention of 
labelling photon energies in the local inertial frame of reference 
by u> , when a connection to quantum processes is most salient, 
and photon energies that an observer at infinity would measure 
will be denoted by e . By the same token, B will be used to 
denote local fields, and B will represent the surface field that 
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would be inferred at the pole in flat spacetime, i.e. a neutron 
star of radius R with a pure dipole field has a magnetic moment 
of m - B Q R 3 /2 . 

Photons of _L polarization dominate the photon population 
generated by the principal primary and secondary emission pro- 
cesses, namely cyclotron/synchrotron radiation, curvature emis- 
sion and resonant Compton scattering. Synchrotron and cur- 
vature radiation are central to the hard X-ray and gamma-ray 
emission of polar cap models (e.g. Daugherty and Harding 1982; 
1996), and are essentially identical in their polarization proper- 
ties in the classical picture. For monoenergetic electrons, polar- 
ization levels V = \h± — n|||/|n_i_ + h\\\ of between 50% and 
100% are achieved (e.g. see Fig. 6.7 of Bekefi 1966), while for 
power-law electrons, the degree of polarization is generally in 
the 60%-70% range (e.g. Bekefi 1966; Rybicki and Lightman 
1979), with _L photons dominating in both cases. At B ;> B CI , 
there is some degree of quantum depolarization, though in re- 
ality this is generally small since photon angles with respect 
to local fields are always small. The dominance of _L photons 
applies also to the products of resonant Compton scattering, 
which has more recently been considered (e.g. Sturner and Der- 
mer 1994) as a primary emission mechanism. In the Thomson 
limit, the cross-sections for scatterings of polarized photons can 
be derived from the results of Herold (1979), indicating that _L 
photons are produced at 3 times the rate of || photons. 



2.1. Competition with Magnetic Photon Splitting 
2.1.1. Photon Splitting 

The relevance of photon splitting 7 — ► 77 to neutron star 
environments was emphasized by Adler (1971), Mitrofanov et 
al. (1986) and Baring (1988). In the context of pulsars, it has 
been discussed by Baring (1993), who proposed it as a mecha- 
nism for suppressing the appearance of e+e~ annihilation lines, 
and by Harding, Baring and Gonthier (1996, 1997) and Chang, 
Chen and Ho (1996), who focused on its action in attenuating 
gamma-ray pulsar continua. Splitting is a third-order QED pro- 
cess with a triangular Feynman diagram. Though it is permitted 
by energy and momentum conservation, when B = it is for- 
bidden by a charge conjugation symmetry of QED known as 
Furry's theorem (e.g. Jauch and Rohrlich 1980), which states 
that ring diagrams that have an odd number of vertices with 
only external photon lines generate interaction matrix elements 
that are identically zero. This symmetry, which pertains to the 
electron/positron propagators, is broken by the presence of an 
external field. The splitting of photons is therefore a purely 
quantum effect, and has appreciable reaction rates only when 
the magnetic field is at least a significant fraction of the quan- 
tum critical field B CI . 

It is practical to restrict considerations of splitting to regimes 
of weak dispersion, where manageable expressions for its rates 
are obtainable, but are still complicated triple integrations (e.g. 
see Adler 1971; Stoncham 1979; Baier, Mil'shtein, & Shaisul- 
tanov 1996; Adler & Schubert 1996) or triple summations (Bar- 
ing 2000). Further specialization to either low magnetic fields 
( B <C B CI ) or low photon energies ( lo <C 2 ) therefore proves 
expedient, and palatable results for splitting rates were first 
obtained in such regimes by Bialynicka-Birula and Bialynicki- 
Birula (1970), Adler et al. (1970) and Adler (1971). A compact 
presentation of these rates (i.e. for lo <C 2 ) for the three po- 
larization modes of splitting permitted by CP (charge-parity) 
invariancc in QED, namely _L— >||||, _L^_L± and ||— >_L|| , is 



(e.g. see HBG97, who also display differential rates) 
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in the frame where photons propagate perpendicular to the field, 
where 
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and Qf is the fine structure constant and A = h/(m e c) is the 
Compton wavelength of the electron divided by 2tt . At low 
fields, Mi and M2 are independent of B, but at high fields 
possess Mi oc B~ 3 M2 oc B~ A dependences. These are the ex- 
pressions used in this paper because of their broad applicability 
to the pulsar problem. Deviations from this low energy limit 
near pair creation threshold are presented in detail by Baier, 
Mil'shtein, & Shaisultanov (1996) and Baring & Harding (1997), 
and are mentioned below where appropriate. An analytic ap- 
proximation to the total rate in the B ;§> B CI limit is presented 
in Baring and Harding (1997), and B ^> B cr analytic forms for 
other splitting modes can be found in Baring (2000). 

The birefringence of the magnetized vacuum implies an al- 
teration of the kinematics of strong field QED processes (Adler 
1971), admitting the possibility of non-collinear photon split- 
ting. Hence, while the splitting modes _l_ — > _l_ 1 1 , ||— >_L_L and 
|| — >|| || are forbidden by CP invariance in the limit of zero disper- 
sion, dispersive effects guarantee a small but non-zero probabil- 
ity for the J.— >-L|| channel. Extensive discussions of linear dis- 
persion in a magnetized vacuum are presented by Adler (1971) 
and Shabad (1975); considerations of plasma dispersion (e.g. see 
the analysis of Bulik 1998) are not relevant to classical gamma- 
ray pulsars because of the relatively low densities present in 
the magnetosphere, but may be quite pertinent to soft gamma 
repeaters. Adler (1971) showed that in the limit of weak lin- 
ear vacuum dispersion (roughly delineated by -Bsin^kB Bcr), 
where the refractive indices for the polarization states are very 
close to unity, energy and momentum could be simultaneously 
conserved only for the splitting mode _L— >|| || (of the modes per- 
mitted by CP invariance) below pair production threshold. This 
kinematic selection rule was demonstrated for linear dispersion, 
a regime that applies to subcritical fields. Therefore, it is prob- 
able that only the one mode (_L— >|||| ) of splitting operates in 
normal pulsars. However, this constraint may not hold in su- 
percritical fields where strong vacuum dispersion arises, thereby 
requiring a revised assessment using the generalized vacuum po- 
larizability tensor (i.e. including quadratic and higher order 
contributions to the vacuum polarization). 
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2.1.2. Attenuation Lengths 

As pair creation is a first order QED process, whereas split- 
ting is third-order in the fine structure constant a.f — e 2 /(he) , 
it is not immediately obvious that 7 — ► 77 can ever domi- 
nate 7 — > e + e~ in pulsar environs. Yet it is the propagation 
of the photons at small angles to the field through the pulsar 
magnetosphere that affords photon splitting an opportunity to 
compete effectively with pair creation when B J> B cr , since the 
pair threshold u> sin 6*kB = 2 is always crossed from below. An 
approximate assessment of the relative importance of magnetic 
photon splitting and pair creation 7 — > e + e~ was performed 
by Harding, Baring and Gonthier (1997) by computing atten- 
uation lengths L for each of these processes. These are the 
scalelengths for attenuation in the neutron star magnetosphere, 
and are defined to be path lengths over which the optical depth 



-(9,e, l)= [ T(9 kB , 
Jo 



u)ds 



(3) 



is unity, i.e. r(0,£, L) = 1 . In computing such attenuation 
lengths, it is essential to fully include the effects of general rela- 
tivity. The reason for this is that the quantum transition rates 
for both splitting and pair creation are strong functions of the 
photon energy and angle and the magnitude of the magnetic 
field strength in the local inertial frame, all of which are strongly 
influenced by curved spacetime. Such an analysis is confined to 
the Schwarzschild metric because the dynamical timescales for 
gamma-ray pulsars are considerably shorter than their period 
(e.g. P = 0.15 sec. for PSR1509-58), so that rotation effects in 
the Kerr metric, for example those due to frame-dragging, can 
be neglected in our photon attenuation analysis. Throughout 
this paper, we assume a neutron star mass, M — 1.4 Mq and 
radius, R = 10 6 cm. 

Values for L are computed assuming that test photons are 
emitted on or above the neutron star surface at some polar co- 
latitude (usually chosen to be at the polar cap rim) and 
propagate outward, initially at a specified angle #kB,o to the 
gravitationally-modified dipole magnetic field; a depiction of 
the geometry is presented in HBG97. Frequently in this pa- 
per, a surface origin of the photons is chosen to provide a con- 
cise and representative presentation of the attenuation prop- 
erties of the pulsar magnetosphere. Such attenuation lengths 
possess power-law behavior at high energies, with L oc e~ 5 / 7 
for photon splitting, and L ~ 2/e for pair creation just above 
threshold (HBG97; proportionalities that hold in both curved 
and flat spacetime). Furthermore, they display the property 
that the attenuation length declines with eolatitude of emis- 
sion, a consequence of the associated increase in curvature of 
the field lines. At low energies, the attenuation lengths diverge 
and photons escape the magnetosphere without attenuation. 
The L — > 00 asymptotes define escape energies (HBG97), be- 
low (above) which spectral transparency (opacity) is effectively 
guaranteed. The existence of such escape energies is a conse- 
quence of the r~ 3 decay of the dipole field: there are always 
sufficiently low photon energies for which the field declines be- 
fore the photons have had sufficient time to attenuate in their 
passage through the magnetosphere. 

2.1.3. Escape Energies 

For each photon trajectory through the magnetosphere, cor- 
responding to a specific set of values for the eolatitude and 
$kB,o of emission at the surface, and for a particular dipole 



surface field Bq , the escape energy e csc is uniquely defined 
and cleanly delineates the energy ranges of source opacity and 
transparency. The escape energy also depends on the radius of 
emission Rq (see Figure ||). The dependence of e csc on such 
parameters is shown in Figure |l| for photons that are initially 
of polarization _L . Escape energy plots for unpolarized photons 
are given in Figures 3 and 5 of HBG97. The feature of Fig- 
ure p] that is most salient for the considerations of this paper 
can be obtained by comparing the plots for splitting and pair 
production. For low fields, pair production escape energies are 
below those for splitting, but the situation is reversed in high 
fields; the escape energies are roughly equal for a narrow band 
of fields around Bq ~ 0.5i? cr . For field strengths Bq f> B cl - , 
photon splitting can attenuate photons well below pair thresh- 
old at significant colatitudes. Hence splitting can be expected 
to dominate 7 — ► e ± in supercritical fields for the attenuation 
of photons of _L polarization. 

The escape energies in Figure [j] generally decline with and 
are monotonically decreasing functions of Bq for the range of 
fields shown. Consider first the dashed curves in each panel, 
corresponding to initial propagation of photons along the field. 
The divergences as — > are due to the divergence of the 
field line radius of curvature at the poles. The maximum an- 
gle 6>kB achieved before the field falls off and inhibits attenua- 
tion is proportional to the eolatitude . For photon splitting, 
since the rate in equation ([!]), and therefore also the inverse of 
the attenuation length L , is proportional to w 5 sin 6 #kB , it fol- 
lows that the escape energy scales as e csc oc -6 / 5 near the 
poles. At fields Bo ;> 0.3-B cr , there is a diminishing depen- 
dence of B in the attenuation coefficient (e.g. see Adler 1971; 
Baring & Harding 1997 for graphical illustrations) so that a sat- 
uration of the photon splitting attenuation lengths and escape 
energies arises in highly supercritical fields. For pair produc- 
tion, the behaviour of the rate (and therefore l/L ) is exponen- 
tial in l/(cjf?sin#kB) (e.g. see Daugherty & Harding 1983), 
which then quickly yields a dependence e csc oc _1 near the 
poles for Bq <J 0.1i? cr . This behavior extends to higher surface 
fields because production then occurs at threshold, which deter- 
mines e osc ~ 2/#kB oc _1 . The pair production escape energy 
curves are bounded below by the pair threshold 2 / sin #kB at 
the point of pair creation (not photon emission) , and for high 
approach the pair threshold (1 + yT+ 2B/B CI )/ sin^kB (with 
sin#kB ~ 1 determined by geometry), blueshifted by the fac- 
tor (1 - 2GM/RC 2 )- 1 / 2 ~ 1.3. The field dependence in this 
saturation energy arises because the creation of pairs by _L - 
polarized photons cannot leave both the electron and positron 
in the lowest Landau level (e.g. see Daugherty & Harding 1983); 
the minimum energy configuration requires one member of the 
pair to be in the first excited state. 

Since pulsar cascade high energy emission from curvature ra- 
diation, inverse Compton or synchrotron by relativistic particles 
with Lorentz factor 7 e will not beam the photons precisely along 
the magnetic field, but within some angle ~ l/7e to the field, it 
is important to illustrate the effect on the escape energies of a 
non-zero angle of emission #kB.o of the photons relative to B. 
This is done via the solid curves in Figure l], where #kB.o = 10~ 2 
(i.e. 0.57° ) is taken towards the dipole axis; qualitatively and 
quantitatively similar behavior arises when #kB,o = 10~ 2 is 
taken away from this axis. Curvature radiation-initiated cas- 
cades generally have 7 e ~ 10 7 (e.g. Daugherty & Harding 1989; 
see also Harding & Muslimov 1998), while inverse-Compton 
seeded pair cascades yield 7 e ~ 3 x 10 5 -10 6 (e.g. Sturner 




Fig. 1. — The escape energy (i.e. where L — > oo ) for (a) photon splitting and (b) pair production as a function of magnetic colatitude for photon 
emission both along B (dashed curves) and at angle $kB,o = 0-01 radians (= 0.57° ) to the field (solid curves). The escape energies for each process are 
monotonically decreasing functions of B for the range of parameters shown. The #kB,o = curves have slopes of -6/5 (splitting) and -1 (pair creation) 
at small ©, as discussed by Harding, Baring & Gonthier (1997), and diverge near 6 = 0, where the field line radius of curvature becomes infinite. For 
the solid curves, at low magnetic colatitudes 6 < 109kB,o i_* ne ne ld curvature is so low that photon attenuation is insensitive to the value of 6 and is 
well described by the uniform field results in Eqs. (hi) and (H). The escape energies are presented for surface photon emission, Rq = R. 



1995; see also Harding & Muslimov 1998), both yielding angles 
0kB,o very much smaller than the example in the figure. The 
situation is similar for synchrotron radiation from first genera- 
tion secondary electrons, where j e ~ 10 3 — 10 4 ; such Lorentz 
factors result from pairs created by primary photons (curvature 
or resonant Compton) in the 1 GeV-10 GeV range. However, as 
the cascade proceeds, higher generations of pairs achieve lower 
Lorentz factors, resulting in a cumulative power-law spectrum 
between -y e ~ 10 2 and j e ~ 10 4 (Daugherty & Harding 1982). 
This behaviour applies to pulsars with low to moderate magnetic 
fields, Bo <J 0.2B CI , a domain that is well-studied in the litera- 
ture. For more highly-magnetized pulsars, the restrictio n of pair 
creation to a single state, discussed at length in Section |2~2] , will 
inhibit the creation of second and higher generation pairs, im- 
plying that the dominant synchrotron signal will sample angles 
to the field much smaller than #kB,o = 10~ 2 . Hence, in summa- 
tion, we expect that #kB,o < 10 -4 — 10 -3 will be representative 
for the considerations of this paper. The #kB,o = 10~ 2 choice 
in Figure [l] aids clarity of illustration. 

The general behaviour of the 6*kB.o = 10~ 2 curves in Figure [l] 
can be simply understood. For the most part, the escape energy 
is insensitive to the emission angle for > 10#kB,o ■ For small 
angles, the escape energy decreases and the curves flatten below 
the #kB,o = curves, converging as 6 — > to an energy that is 
proportional to (B sin#kB,o)~ 6 / 5 when B a <C B cr . This con- 
vergence is a consequence of the field being almost uniform and 
tilted at about angle #kB,o to the photon path for trajectories 
that originate near the pole; this behaviour at low colatitudes 
was first noted, in the case of pair creation in flat spacetime, by 
Chang, Chen and Ho (1996). These "saturation" energies fol- 
low from the dependence of the photon splitting rate on B and 
, and have a declining sensitivity to Bq when Bo ;> 0.3B cr . 
In Fig. gb, the same effect is seen for pair creation, but this 
time the "saturation" is at the rcdshiftcd threshold energy for 
7 — > , which varies as (1 + yl + 2£?/i? cr )/ sin#kB,o • Hence, 



instead of the pair creation escape energies always being mono- 
tonically decreasing functions of Bq , they experience an inver- 
sion at small colatitudes and increase with Bq . Such an effect 
was not present in the polarization-averaged escape energies pre- 
sented in HBG97 because the || state could always access the 
true pair threshold, i.e. 2/ sin^kB.o ■ Note also, that while split- 
ting escape energies always drop when #kB,o is increased from 
zero to 10~ 2 , opposite behaviour is seen in high fields for pair 
creation at moderate colatitudes, due to subtleties concerning 
the sudden onset of pair creation (HBG97). 

It is important to emphasize that general relativistic effects 
are crucial to these calculations. While the dependence of the 
escape energy on emission colatitude and surface polar field 
strength is similar in curved and flat spacetime, the introduction 
of the Schwarzschild metric results in decreases of the escape en- 
ergies for both process by factors of between 2 and 4, as is illus- 
trated in Figure 4 of HBG97. Moreover, these decreases differ 
for splitting and pair creation, making it imperative to carefully 
account for propagation in curved spacetime when comparing 
attenuation characteristics for the two processes. The largest 
effects of the incorporation of general relativity are due to the 
increase of the surface dipole field strength by roughly a fac- 
tor of 1.4 above the dipole spin-down estimate of Bo , and the 
correction for the gravitational redshift of the photon, which 
increases the photon energy by roughly a factor 1.2-1.3 in the 
local inertial frame at the neutron star surface compared to the 
energy measured by the observer in flat spacetime at infinity. 
The influence of these modification factors is amplified by the 
sensitivity of the rates for 7—5-77 and 7 — > e + e~ to B and e 
so that 20%-30% deviations from flat spacetime parameters map 
over to the factors of 2-4 in the escape energies. Near the polar 
cap, the curvature of the photon trajectory in a Schwarzschild 
metric does not affect the escape energies, to first order, since it 
is generally compensated for by relativistic modifications to the 
curvature of the magnetic field (e.g. see Gonthier & Harding 
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Fig. 2. — The solutions to Eq. M) that define (a) the polar surface field Bo and (b) the escape energy e esc (for either photon splittings _L— >|||| or pair 
creation _L— > e + e~ ) as functions of the colatitude © of the point of photon emission, which is assumed to take place on the stellar surface. Solutions are 
displayed for the four values of S^b o i as labelled, with the $kB,o = case generating the approximate dependences Bo oc -4 / 15 and e eS c oc _1 . All 
$kB,o 7^ cases produce B and e e sc that are independent of © for sufficiently small colatitudes, namely roughly when < 10 4 9] i b,o ■ 



1994). 

The most concise designation of when the creation of pairs 
might be suppressed by photon splitting is when the escape en- 
ergies of the two processes are set equal: 



sir 1111 ^, e, Ro) = 4r e+e ~(£o, e, R ) 



(4) 



where Rq is the radius of the point of emission. Note that 
defining a similar equality for the || polarization mode, or a 
polarization-averaged alternative, leads to only very modest 
changes in the ensuing results. Baring & Harding (1998) con- 
tended that Eq. (H) is representative of the delineation between 
pulsar parameter regimes where pairs are created in profusion 
and when they are present only in paucity (i.e. at higher Bq ). 
The purpose of the extensive investigation of this paper is to 
determine under what conditions such a demarcation is truly 
representative of a suppression of pair creation in pulsars; such 
an exploration appears in Section 3.2 below. The solutions of 
this equality define curves in the (B , 0) plane, and therefore 
specify functional dependences of the dipole spin-down mag- 
netic field Bq and the escape energy (for either process) on the 
colatitude of emission. These functions are displayed in Fig- 
ure H for different values of the initial angle #kB,o of the photons 
relative to the field. These solutions define the basis for our con- 
siderations of the suppression of pair creation in the context of 
radio pulsars in Section || below. 

To understand the behaviour of the curves, first consider the 
#kB,o = case. The escape energies plotted in Fig. [j] have the 
approximate dependence e osc oc B^ a , with a w 3/4 (for 
B cl < Bq < &B CT ) and (3 = 6/5 for splitting, and (3=1 for 
pair creation. The dependence of the pair creation escape energy 
on field strength for Bq ;> B cl - is quite weak ( a « ) so that 
solutions to Eq. (Q) must yield escape energies approximately 
proportional to _1 . This can be immediately folded into the 
photon splitting proportionality e csc oc Bq 3 ^ 4 0~ 6 / 5 to yield 
the magnetic field dependence Bq oc 0~ 4 / 15 in Figure^, an ap- 
proximation that is applicable for B cr <^ Bq < 4 B cr , the range 



of fields of interest to the focus of this paper. The variation of 
Bq with actually increases very slowly at smaller colatitudes 
due to the diminishing dependence of the photon splitting rate 
for the J-— >|||| mode on field strength in highly supercritical 
fields. The non-zero #kB,o solutions for Bq are monotonically 
decreasing functions of 0kB.o i an d eventually become indepen- 
dent of colatitude when <; 10 4 #kB,o , as is expected from the 
horizontal branches of the escape energy plots in Figure |2|. Since 
the escape energy solutions trace the behaviour of the pair cre- 
ation e csc , small increases with #kB,o are exhibited at moderate 
to high colatitudes, even though for <; 10 4 6*kB,o j the solutions 
scale as (6*kB,o) _1 , as expected. 

2.2. Threshold Pair Production 

As discussed in the previous sections, pair production by pho- 
tons initially emitted at small angles to the field occurs very 
near threshold in magnetic fields exceeding 0.1i? cr . In this case, 
there are only a small number of kinematically available electron 
and positron Landau states. In very high fields, the accessible 
number of excited pair states diminishes, and the pairs created 
by most photons in the primary particle spectrum will occupy 
the ground state (for photon polarization || ) or the first excited 
state (for photon polarization _L ). Pairs in the ground state can- 
not produce synchrotron photons. Even pairs in the first excited 
state radiate photons of energy u ~ + 2B/B cr + pi — 1 for 
electron momentum p z parallel to the field, which will gener- 
ate very few pairs (Harding & Daugherty 1983), and certainly 
none if p 2 z < 4(1 + yT-j- 2B/B CI ) for photons of _L polariza- 
tion. This phenomenon thus provides another mechanism for 
suppressing the number of pairs produced by pulsar cascades. 
Although pairs in the ground state could be excited to higher 
Landau levels through inverse-Compton scattering of soft X-ray 
photons (e.g. Zhang & Harding 2000), we leave discussion of 
the potential importance of this effect to the end of this section. 

First, we examine the effect of threshold pair production. 
The threshold photon energy in the center-of-momentum (CM) 
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frame for excitation into the electron-positron Landau state 
combination (j,k) is 



e hk = y/i + 2jB/B CI + y/i + 2kB/B c 



(5) 



As was discussed at the beginning of this Section, the photons 
radiated in pulsar cascades are predominantly of _L polarization, 
which have a pair production threshold at (0,1), i.e. one mem- 
ber of the pair is in the ground state and the other is in the first 
excited state. In a uniform magnetic field B , the number of 
pair states available to a photon which pair produces at energy 
lu cm = uj sin 6>kB in the CM frame (for propagation angles #kB 
to the field), is approximately (Daugherty & Harding 1983) 



1 f B \ ^ <■ 

-^states (w CM , B) Ri ^4 ( J W CM (w OM + 4) (w CM — 2) z 



w sin fkB 



We can assess the importance of ground-state pair creation by 
computing N states (uj suiO^q, B) at the point of pair cr eation 



in the attenuation length calculation described in Section 2.1.2 



Photons were injected at the stellar surface with momenta along 
the local field lines and followed along their curved trajecto- 
ries through the magnetosphere. For very low photon ener- 
gies uj sin 9 kB < 1 + (1 + 2j max B/B„) 1 ^ 2 , A sta tc S (w sin6> kB , B) 
is sampled from a table constructed by summing the exact 
number of pair states (i.e. threshold energy values ujj^ ex- 
ceeded) as a function of cjsin^kB • At higher photon energies, 
-States sin (9kB, B) is computed directly from equation (^) for 
the particular values of uj sin t?kB that are sampled at the point 
of pair production. 
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Fig. 3. — The number of kinematically-accessible electron-positron pair 
states Nutates , at the point of pair creation, of a photon having polarization 
state _L , emitted from the neutron star surface, parallel to the magnetic 
field at colatitude © = 5° . Abates is displayed as a function of the photon 
energy in the observer's frame at infinity. The curves are labeled with val- 
ues of the dipole spin-down magnetic field strength Bo . Single-state pair 
creation becomes rife when Bq > 0.15-B C r ■ The vertical asymptotes to the 
left mark the respective escape energies. 

Figure || shows the number of available pair states at the 
point of pair production as a function of the lab-frame energy 
of a photon emitted parallel to the magnetic field at colatitude 



9 = 5°. Essentially, this Figure is a modified version of Figure 2 
of Harding & Daugherty (1983), which treated only flat space- 
time and did not examine behaviour near the escape energies 
(i.e. regimes where field decline with radius plays a role). The 
curves in Figure || all increase sharply at low photon energies, 
near the escape energy, where the attenuation length goes to 
infinity. At the escape energy, N statcs (u> CM , B) formally goes 
to infinity because photons at this energy have large attenua- 
tion lengths and thus reach the outer magnetosphere where the 
field strength has decreased. At high energies, pair creation 
accesses many states even though it occurs in the inner magne- 
tosphere, and the limit ajsinc?kB ^ 1 is realized. In such cases, 
the attenuation length is L oc p c #kB > since the path length 
s scales as s « p c @kB for radius of curvature p c (oc G -1 ), 
and the function T^ p (w sin#kB) acts like a step function in the 



(6) spatial integration determining L . Hence #kB oc p t 



-1/2 



and 



the number of accessible pair states asymptotically approaches 
•^states oc e 4 Bq 2 p~ 2 oc e 4 2 / B 2 . Therefore, an increase in the 
field Bo or a decrease in the polar cap colatitude 9 reduces the 
number of accessible pair states. Because the pair attenuation 
coefficient increases exponentially as exp(— AB ci /[3loB sin ^b]) ? 
pair production occurs well above threshold in relatively low 
surface fields ( Bq <C 0.1i? cr ), and photons of all energies create 
pairs in highly excited states. These pairs then produce copious 
numbers of synchrotron photons, which initiate pair cascades 
with multiplicities of ~ 10 3 — 10 4 (e.g. Daugherty & Harding 
1982). In higher fields (Bq > . 1 i? cr ) , some photons in the 
spectrum create pairs in the ground state (0,1) or (1,0) for _L 
polarization. By Bq = 0.5B CT , virtually the entire spectrum of 
photons will create pairs in state (0,1) or (1,0). The "noisiness" 
at the minimum of the Bo =0.1 curve in Figure || is a real dis- 
creteness effect, due to the fact that the field at the pair creation 
point is increasing with photon energy and the ordering of the 
(j,k) states depends on the local field strength. From these re- 
sults, we can conclude that there will be a decrease in pair yield 
in pulsar polar cap cascades at magnetic fields Bo > 0.1B cr . 
For Bo > 0.2B cr , there will be strong suppression of cascade 
pair yields. Note that this condition relates to surface polar 
fields. The site for pair creation is generally proximate to the 
pole only for photon energies exceeding around 1 GeV in Fig- 
ure [| otherwise the local field drops below Bo as the escape 
energy is approached. 

As mentioned above, inverse-Compton scattering of pairs in 
the ground state could excite them to higher Landau states, 
thereby increasing the number of synchrotron photons radiated 
and subsequently the pair yields. In the case of highly rela- 
tivistic particles moving along the magnetic field, the incident 
photons have angles to the field near zero in the particle rest 
frame. For such photon incidence angles, the scattering cross 
section in a magnetic field is strongly suppressed below the cy- 
clotron resonance, lu b = B / B cr . Above the resonance, the par- 
ticle may be excited to a higher Landau state. Although the 
scattering cross section does not have resonances at photon en- 
ergies above the cyclotron fundamental for incidence along the 
field, excitation through continuum scattering at these energies 
is still possible, and becomes more probable with increasing field 
strength (Daugherty & Harding 1986; Gonthier et al. 2000). If 
the incident soft photons have a quasi-isotropic blackbody spec- 
trum at temperature T and their average energy is uj s — 3kT , 
then scattering above the resonance in the particle rest frame 
will occur when 7 £ (B/B CT ) m e c 2 /3kT = 2 x 10 3 (B/B CI ) Tg 1 , 
a condition applicable for all fields provided that the particles 
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are fairly relativistic. The mean- free path of these particles to 
scattering will then be 



X, 



1 



7.5 x 10 4 ?T T 6 cm, 



(7) 



where T 6 = T/10 6 K and 77 is a suppression factor for the 
scattering cross section a = r/ar in the Klein-Nishina regime. 
Such mean free paths are short enough to suggest that exci- 
tation through scattering is potentially important. However, 
these photons scattering above the resonance in the electron rest 
frame will not form the bulk of the produced pairs, because the 
resonant cross section (and thus the scattering rate) is several 
orders of magnitude larger than for non-resonant cross section. 

2.3. Positronium Formation 

Another mode of pair creation exists, namely the formation 
of pairs in a bound state, i.e. positronium. This has been pro- 
posed as an effective competitor to the production of free pairs 
(Shabad & Usov 1985, 1986; Herold, Ruder & Wunner 1985; 
Usov & Shabad 1985; Usov & Melrose 1995) because the bind- 
ing energy lowers the threshold slightly ( 1 %) below the value 
for production of free pairs. Positronium formation fundamen- 
tally alters the magnetosphere: it amounts to the suppression of 
electron or positron currents that can screen the induced electric 
fields in the rotator. Furthermore, depending on how long the 
pairs remain bound in their passage along open field lines to the 
light cylinder, such a presence of neutral positronium will help 
to suppress any collective plasma modes that might spawn radio 
emission. Hence positronium formation is a viable alternative 
means of inhibiting the radio signal from pulsars. 

The relevance of positronium formation to pulsars is poten- 
tially great if the bound state is stable for considerable times. 
Positronium is subject to destruction by three main mecha- 
nisms: free decay, electric field ionization, and photo-ionization. 
Using the spontaneous two-photon decay rate (in the positron- 
ium rest frame) of 3.5 x 10 14 (-B/-B cr )sec _1 computed by Wun- 
ner & Herold (1979), Bhatia, Chopra & Panchapakesan (1987) 
obtained positronium decay lengths considerably shorter (by 
two orders of magnitude for bulk Lorentz factor 7 P = 10 3 ) 
than the stellar radius for B ~ 0.2£? cr ; the length is a strongly 
declining function of B . Hence it appears that positronium is 
unlikely to be long-lived in the magnetosphere except for a nar- 
row range of fields, or for sufficiently high positronium Lorentz 
factors, i.e. from the attenuation of photons above 10 GeV. Ap- 
proximate estimates of when field ionization becomes important 
were obtained by Herold, Ruder & Wunner (1985) and Usov & 
Melrose (1996): for pulsar periods P £ 0.5(10 J B / J B cr ) 2 / 3 sec, 
ionization due to E\\ should convert positronium into free pairs. 
This corresponds to a sizeable portion of the phase space of in- 
terest for the considerations of this paper. In addition, photo- 
ionization can destroy positronium, forming free pairs. Herold, 
Ruder & Wunner (1985) claimed that this was highly likely due 
to collisions with thermal radiation emanating from the stel- 
lar surface for temperatures greater than around 10 5 K. Usov 
& Melrose (1995) used the calculations of Bhatia, Chopra & 
Panchapakesan (1992) to argue that photo-ionization timescales 
exceed dynamical ones for surface temperatures ~ 10 5 K. This 
discrepancy is yet to be resolved. 

2.4. Pair Cascades in High Magnetic Fields 

We have discussed several of the physical mechanisms capa- 
ble of suppressing pair production when the local magnetic fields 



approach and exceed the critical field. To more realistically mea- 
sure pair suppression in high magnetic fields, one must examine 
not only the propagation of single high-energy photons through 
the neutron star magnetosphere (i.e., the escape energies) but 
also the transport of the higher generations of photons and par- 
ticles, which are produced by the primary photons that do not 
escape. Even though the first generation of photons split in- 
stead of producing pairs, the second generation of photons may 
create pairs. It is therefore necessary to investigate pair yields 
of the pair/photon ca scad es which are initiated by the primary 
photons. In Section 3.2, we will present results of numerical 
simulations of such pair cascades. In this section, we will qual- 
itatively discuss how we expect the nature of these cascades to 
change as the magnetic field increases. 



LOW FIELD CASCADES 

B < 0.1 B„ 



No Splittin 



HIGH FIELD CASCADES 

B > 0.1 B„ 



3 Splitting Modes 
1 -> II II 1 -> II II 

Splitting Only 1 -> 11 

II -> 111 




Fig. 4. — The nature of pair cascades in increasing magnetic fields. 
Wavey lines represent photons and straight lines represent electrons and 
positrons. Polarization modes of photons ( || - parallel or _L - perpendicu- 
lar) are noted where relevant. 

Figure ^ illustrates the behavior of low and high field pair cas- 
cades. Pair cascades in magnetic fields well below B CI have been 
extensively studied in connection with models of pulsar high en- 
ergy emission (e.g. Daugherty & Harding 1982, 1996). In fields 
B <^ 0.1-B cr , both _L and || mode photons above escape energy 
will produce pairs well above threshold, so that the resulting 
electrons and positrons will initially occupy high Landau states. 
Each pair member will make multiple Landau state transitions 
before reaching the ground state, radiating many synchrotron 
photons. Many of the synchrotron photons will pair produce, 
initiating more synchrotron emission sequences. Such a process 
can go on for a number of generations, provided that the energy 
of the primary photon is high enough, resulting in the produc- 
tion of many pairs for each primary photon. 

When B > 0.1B CI , the prima ry p hotons create pairs at pair 
threshold, as discussed in Section 2J2 ■ This has a profound effect 
on the pair cascades, since the || mode photons produce pairs 
with both members in the ground state, cutting off all further 
generations. Photons of _L mode, which have a higher thresh- 
old, produce pairs with one member in the first excited state. 
The single Landau state transition made by this particle will 
result in primarily _L mode photons which may create a sec- 
ond generation pair with one member in the first excited state. 
Higher generations of pairs are therefore possible, depending on 
the energy of the primary photon, but the total pair yield is 
severely diminished because there is no multiplication of pairs 
with each generation, as in the low-field cascades. 
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At fields B > 0.5B cr , the photon splitting attenuation 
lengths become shorter than pair attenuation lengths. At this 
point, the nature of the pair cascades will depend on which 
modes of splitting are operating at what field strengths. In 
fields B <J B CI , it is probable that only the _L — >-|| || mode op- 
erates (see Section 2.1.1). Thus, primary || mode photons can 
produce pairs, but with both members in the ground state. The 
primary _L mode photons above photon splitting escape energy 
will split into two || mode photons, each of which can produce 
a pair in the ground state. Thus, each polarization branch of 
the cascade ends after, at most, one pair generation. The action 
of photon splitting in one mode therefore will lower the cascade 
pair yield above 0.5-B cr , but only mod erately. 

If, as was discussed in Section 2.1, other modes of photon 
splitting open up in magnetic fields above B cr due to increased 
vacuum dispersion, then the character of pair cascades undergo 
a further transition. If both photon polarization modes can 
split (as is the case when all three modes permitted by QED arc 
operating), then it is possible to completely prevent the produc- 
tion of pairs. A pure photon splitting cascade then ensues, in 
which photons split repeatedly until they can escape the magne- 
tosphere. The only question then becomes whether the splitting 
cascade can degrade the photon energy fast enough so that the 
photons escape without pair conversions in lower fields. Since 
the photons are splitting below pair threshold in high fields, this 
would seem very likely. 

3. PAIR SUPPRESSION IN HIGHLY-MAGNETIZED PULSARS 

Pair creation is obviously pertinent to the hard X-ray and 
gamma-ray emission of pulsars. Yet it is also relevant to the 
discussion of coherent emission in radio pulsars, since it is com- 
monly assumed that a plentiful supply of pairs is a prerequisite 
for, and maybe also a guarantee of, coherent radio emission at 
observable flux levels. Such a connection is the premise of stan- 
dard models for radio pulsars (e.g. Sturrock 1971; Rudcrman 
& Sutherland 1975; Arons & Scharlemann 1979), a relationship 
founded in the simplicity of its explanation (Sturrock, Baker & 
Turk 1976) for the extinction of radio pulsars beyond the con- 
ventional death line at long periods. Detailed discussions of the 
relationship of pairs to the production of radio emission can 
be found in Michel (1991) and Melrose (1993), though we note 
that there are dissenting views to this popular connection (e.g. 
Weatherall & Eilek 1997). 

The pre-eminent consequence of any suppression of pair cre- 
ation in pulsars is that the emission of radio waves should be 
strongly inhibited. Hence the issue of possible quenching at 
high Bq by the mechanisms discussed above forms the focus of 
this Section. First we compute the putative radio quiescence 
boundaries in the P-P diagram using the approximate cri- 
terion for the sup pressi on of pair creation by photon splitting 



3.1. Radio Quiescence Boundaries in the P-P Diagram 



defined in Section |2.1.3, namely when the escape energy for the 



splitting of _L photons is less than or equal to that for pair cre- 
ation _L— > e^- . Essentially this provides a detailed derivation of 
the boundary introduced in Baring & Harding (1998). We then 
model the physics of pair suppression near pulsar polar caps 
by means of a Monte Carlo simulation of pair cascades in high 
fields, treating both photon polarization states. These detailed 
calculations identify the requirements for effective pair suppres- 
sion which must be met t o just ify the simplified criterion for 
pair suppression of Section 2.1.3. 



In Section 2.1.3, the condition for the equality of photon split- 



ting and pair creation escape energies established a region in 
(Bq, 0) space where pair suppression would probably occur. 
This space can be easily transformed into a region in the ra- 
dio astronomer's traditional pulsar phase space, which consists 
of the measurables P , the pulsar period, and P , the period 
derivative. The first part of this transformation derives from 
the relationship between the period and the size of the polar 
cap . Such a relation follows from the definition of the cap as 
the portion of the stellar surface that anchors open field lines, 
i.e. those that are not closed within the light cylinder. In flat 
spacetime, for a rotating dipole field, the angle 26 subtended by 
the polar cap can be expressed via (e.g. Manchester and Taylor 
1977) sin6 = [2ttR/(P c)} 1/2 . Due to general relativistic dis- 
tortions of the dipole in a Schwarzschild metric, this formula is 
modified to the form (Muslimov and Harding 1997) 



. Q J 2-kR _/ Rs 



1/2 



T{X) = --g" 



-1 -1 



log e (l -*) + §(* + 2) 



(8) 



that is used in this paper. Here R s — 2GM/c 2 is the 
Schwarzschild radius, and the result for any radius r can be 
obtained by the substitution R — > r . Clearly, the strong gravi- 
tational field reduces the size of the polar cap (J-(x) 1 — 3x/4 
for i«l), arising in conjunction with its intensification of the 
magnetic field. Note that strictly speaking, the polar cap con- 
nects to field lines that are not closed within the Alfven radius 
r A , however, the intense pulsar fields yield relativistic Alfven 
speeds so that r A is generally close to the light cylinder ra- 
dius P cj (2ir) . The loading of the field with plasma can have 
a significant impact on r A and in soft gamma repeaters, as 
discussed in Section 4.2 below. 



The second half of the transformation of the (Bq, 0) space 
to the (P, P) phase space arises from the commonly assumed 
spin-down relationship between the dipole spin-down field Bq 
and the measurables. For radio pulsars, the increase in period 
is usually attributed to magnetic dipole radiation that taps the 
rotational kinetic energy of the neutron star. This couples Bq 
to P and P through equating the rotational energy loss to 
the dipole radiation energy loss which is specified in terms of 
the total magnetic moment /i B • The radio pulsar community 
has commonly adopted the approximation fj, B ~ BqR 3 (e.g. 
see Manchester and Taylor 1977). The actual dipole moment 
of a star-centered dipole, regardless of internal field configu- 
ration, is fi B — (BoR 3 )/2 (e.g. see Shapiro and Teukolsky 
1983, Usov and Melrose 1995). This choice, which we adopt 
throughout this paper, leads to a dipole radiation loss rate of 
dE/dt — — Bq i? 6 57 4 /(6c 3 ) for an orthogonal rotator of angular 
frequency SI = 2ir/P , and 



B Q = 6.4 x 10 1 



Gauss 



(9) 



for the surface field at the magnetic pole. This estimate is twice 
the conventional choice of Manchester and Taylor (1977). We 
note that multipole contributions near the stellar surface and the 
field geometry in the outer magnetosphere can also influence the 
pulsar surface fields derived from P and P . The uncertainties 
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in the understanding, of the global field structure are sufficient 
that the use of Eq. (g) is justified for the purposes of this paper. 
Note that curved spacetime does not impact the determination 
of B Q , a flat spacetime quantity, but just increases the value of 
the field in the local frame at the pole over Bq ; this increase in 
magnetic field energy density can be viewed as a general rela- 
tivistic "redshifting" of the Poynting flux of the rotator. 
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High B Radio Quiescence Boundaries 
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Fig. 5. — The conventional depiction of pulsar phase space, the P-P 
diagram, with filled circles denoting the locations of 541 members (those 
with P ^ n) nf thr. latcat crlHinn (Taylnr Manchester & Lyne 1995; see 
also http://pulsar.princeton.edu/) of the Princeton Pulsar Catalogue, 



and open circles marking 199 pulsars in the recent Parkes Multi-Ream sur- 
vey irttp : //www. atnf . csiro . au/~pulsar/psr/pmsurv/pmgww/[. Three of 



the handful oi gamma-ray pulsars, namely the tjrab, Vela and PSR 1509- 
58 are highlighted, with point styles as indicated in the inset. The dotted 
diagonal lines denote constant field strength, as inferred from Eq. (H) 
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icial positioning nf the conventional death line, from P-P diagrams 
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it. The heavy solid curves 



is depicted as the dashed line on the 
lve a variety of choices for the boundary of 
radio quiescence for surface emission, depending on the value of SjtB > the 
initial angle of photons relative to the field. In addition, the heavy dashed 
curve is a similar boundary for #kB.o = an d emission point half a stellar 
radius above the surface (i.e. Ro = 1.5R). Each of these curves, defined 
by solutions to Eqs. M)— (H), purportedly partitions the P-P diagram into 
regions of radio loud^below the curve) and radio quiet (above) pulsars, 
the latter being where photon splitting suppresses pair creation for the _l_ 
polarization state. Values 0kB,o < 10 -4 and r > 1.2R can comfortably 
accommodate the current radio pulsar population. 

By inverting Eq. (|J) to solve for P in terms of B and P , 
the solutions Bq(Q) to Eq. (Eh that demarcate the domain of 



possibly strong suppression of pair creation by photon splitting 
are mapped onto the P-P diagram in Figure |^. Curves are de- 
picted for three of the four values of #kB,o that are addressed in 
Fig. ^, specifically for the case of surface photon emission, and 
also for a #kB.o = case of emission above the surface. Much 
of the extant isolated radio pulsar population is also exhibited 
(only those with positive P), namely those sources listed in 
the Princeton Pulsar Catal ogue (Taylor, Manchester fc Lyri c 
1993, with 541 pulsars: see frittp : //pulsar . princeton . edu7 ), 
with a recent update from the Parkes Multi-Beam Survey 
(adding 122 new pulsars with presently archived P and P : s ee 



ittp : //www. atnf . csiro . au/~pulsar/psr/pmsurv/pmwww/ ). 

Note that, from the discussion in Section |2.1.3| , the 6>kB,o 
boundary for surface emission ( Rq = R) corresponds roughly 
to Bq oc O -4 / 15 and hence to F a p-n/15 . specifically, we 
find the approximation 



P « 7.9 x 1(T 13 



P 
1 sec 



-11/15 



(10) 



the slope of which depends on the how strongly the rate of pho- 
ton splitting depends on B in this regime of fields. It should also 
be noted that while such putative boundaries of radio quiescence 
have been labelled by specific values of 0kB,o for the purpose 
of illustration, in fact the value of this initial angle of photon 
propagation with respect to the field is a weak function of the 
colatitude. This is due to the coupling between the Lorentz fac- 
tor of accelerated electrons (and therefore 6>kB,o ) and the polar 
cap size , a connection discussed below in Section |4.2| . Also 
depicted in Figure || is a fiducial positioning (adopted by Tay- 
lor, Manchester & Lyne 1993) of the conventional death line, 
with P oc P 3 (i.e. corresponding fixed open field line voltage, 
adopt), to the right of which there is just one detected radio 
pulsar (PSR J2144-3933). It is a well-known problem that the 
computed position of the pulsar "death line" , assuming pairs are 
produced only by curvature radiation photons, lies at smaller 
periods than the observed "death line" , for the assumption of 
the standard polar cap size (e.g. Arons & Scharlemann 1979). 
However, increasing the polar cap size to less than twice that of 
the standard brings the computed "death line" into agreement 
with virtually the entire radio pulsar population, the notable 
exception being the recently "re-discovered" 8.5 second Parkes 
pulsar (Young, Manchester & Johnston 1999) that is conspicu- 
ous on the right hand side of Fig. ||. However, pair production by 
inverse Compton-scattered photons can explain radio emission 
from this pulsar (Zhang et al. 2000). 

The Princeton catalog by itself clearly rules out boundaries 
of radio quiescence due to photon splitting that correspond 
to #kB,o <; 3 x 10~ 4 for surface emission. The fact that the 
0kB,o = , Rq = R boundary was comfortably located above 
the entire Princeton collection of radio pulsars was an attractive 
feature that underpinned the suggestion of Baring & Harding 
(1998) that this particular case marked the approximate loca- 
tion of the boundary of radio quiescence, where photon split- 
ting could strongly inhibit pair creation. However, this conclu- 
sion has been challenged by the dramatic increase in the ob- 
served population due to the exciting new Parkes Multi-Beam 
survey (Camilo et al. 2000; D'Amico et al. 2000; Kaspi et 
al. 2000), with three new pulsars (PSRs J1119-6127, 1726-3530 
and 1814-1744, all with high dispersion measures in the 700-850 
range) having been discovered that are obviously at P above 
the 6*kB,o = , Rq = R boundary. 
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Although the boundary of radio quiescence for surface emis- 
sion (chosen for illustrative purposes by Baring & Harding 1998) 
neatly falls between the population in the Princeton catalog and 
the handful of known anomalous X-ray pulsars and soft gamma 
repeaters (see the blow-up of the P-P diagram in Figure ^|), 
most (or perhaps all) of which are radio quiet, clearly the new 
Parkes survey data questions this assumption of surface emis- 
sion. The quiescence boundary rapidly moves up on the P- 
P diagram for emission at higher altitudes, which might nat- 
urally be expected for traditional curvature radiation-initiated 
cascades in pulsars with Crab-like surface fields (e.g. around R 
above the surface for Daugherty & Harding's 1996 modeling of 
the Vela light curve). However, for the high B fields sampled 
by the #kB,o <S 10~ 4 curves, inverse Compton is probably more 
relevant to the acceleration region (see Sturner 1995; Harding & 
Muslimov 1998) and low alt itud e emission is much more prob- 
able, as discussed in Section 4.2 below. One might then expect 
that the mean radius of emission satisfies r <J 1.5R. Accord- 
ingly, we have computed a #kB,o = 0, Rq = 1.5R boundary 
and depicted it in Figure ||; it clearly lies above even the highly- 
magnetized Parkes Multi-Beam survey pulsars. Intuitively, it 
might be anticipated that the location of the boundary as a 
function of emission radius Ro should scale simply by the rela- 
tionship between field strength and radius for a dipolar struc- 
ture, so that R>o oc r 3 . This appears to be borne out near the 
stellar surface, however the radial dependence of the location of 
the quiescence boundary weakens at altitudes r ~ R (in spite 
of the reduced effects of general relativity at r > R), princi- 
pally because of changes in the flaring of field lines sampled by 
photons in their flight away from the neutron star. The present 
Parkes survey data constrain computed quiescence boundaries 
to the r ;> 1.2R range; in the light of the present indetermi- 
nacy in th e ac tual locale of the acceleration region (addressed 
in Section 4.2), this is not a serious problem. 



3.2. Simulations of Photon Splitting/Pair Cascades 

Using the equality of photon splitting and pair creation escape 
energies for _L polarization as the criterion for when splitting 
starts to strongly suppress the production of pairs is a simple 
choice, motivated largely by the predominance of production of 
_L photons in the relevant emission processes. It is necessary to 
determine whether and under what conditions this choice is an 
appropriate one. 

To assess the applicability of the radio quiescence boundaries, 
we have studied pair suppression in high magnetic fields by 
means of Monte Carlo simulations of photon splitting/pair cas- 
cades in a neutron star magnetosphere. This calculation gener- 
alizes that of HBG97 by including the cyclotron and synchrotron 
radiation from the electron-positron pairs. This improvement is 
necessary in order to compute the pair yields from all gener- 
ations of the cascade. Note that we omit resonant Compton 
scattering from our simulations since, while it is important as 
a radiation process for primary and secondary particles, it will 
not significantly impact the spatial transfer of photons in normal 
radio pulsars due to the low densities and small angles; it may, 
however, be significant in soft gamma repeaters. We inject pho- 
tons parallel to the local magnetic fields at the neutron star sur- 
face with specified magnetic colatitude and four-momentum 
k . The photon energies are sampled from a power-law distribu- 
tion, 

N(e) = N E~ a , e min < e < e max , (11) 
where we generally set e m i n = 1 , and let e ma x be a free param- 



eter. In the results presented in this paper, we generally take 
a = 1.6 as a spectral index representative of either primary cur- 
vature radiation with energy losses (giving a = 5/3) or resonant 
inverse Compton emission which is relatively hard. Polarization 
is chosen randomly such that the average distribution has 75% 
in the _L mode and 25% in the || mode, reflecting the po- 
larization produced by the radiation mechanisms discussed at 
the beginning of Section ||. The path of each input photon is 
traced through the magnetic field, in curved spacetime, accu- 



mulating the survival probabilities for splitting, P s s , 
pair production, P[ urv , independently: 



Psurv(s) = expj-T(8,£, l)\ 



and for 



(12) 



where r(9,e, I) is the optical depth as defined in Eqn. (||). 
Each photon may split, produce pairs or escape, based on a 
combination of the running survival probabilities for splitting 
and pair production (see HBG97). If a photon splits, the ener- 
gies and polarizations of the final photons are sampled from the 
distribution from the branching ratios given in HBG97. Each 
final photon is then followed in the same way as the parent 
photon. If a photon produces pairs, the total energy, Landau 
state and parallel momentum of the electron and positron are 
determined. Each member of the pair is assumed to have half 
the energy and the same direction of the parent photon (tech- 
nically this is only appropriate for high Landau states), except 
when the pair is produced at threshold (i.e. in the [0,0] state for 
|| polarization and [0,1] or [1,0] for _L polarization), in which 
case we use the full kinematic equations to determine the en- 
ergy and momentum of each pair member. Each member of the 
pair occupying an excited state emits a sequence of cyclotron 
or synchrotron photons. The method used to simulate the cy- 
clotron/synchrotron emission is similar to that of Daugherty & 
Harding (1996). If the particle Landau level is larger than 20, 
the high-energy limit of the quantum synchrotron transition rate 
(Sokolov & Ternov 1968) is used, in which case we assume that 
the photons are emitted perpendicular to the magnetic field in 
the particle rest frame (high-energy limit). When the particle 
Landau level is smaller than 20, the exact QED cyclotron tran- 
sition rate (Harding & Preece 1987) is used, in which case the 
angles of the emitted photons are sampled from a distribution. 
In both cases, the emitted photon polarizations are sampled 
from the corresponding polarization distributions. Each emit- 
ted photon is propagated through the magnetic field from its 
emission point until it splits, produces pairs or escapes. The 
cyclotron/synchrotron emission sequence continues until each 
particle reaches the ground state. The cascade continues until 
all photons from each branch have escaped. Throughout the 
cascade, a running tally is kept of the number of pairs produced 
and of the number of photon splittings. General relativistic ef- 
fects of a Schwarzschild metric on the photon momentum and 
dipole magnetic field, are included in the same way as described 
in HBG97. 

We have chosen the total number of pairs produced per in- 
jected photon as a quantitative measure of the pair yield of the 
cascade initiated by a particular parent photon spectrum. The 
free parameters are then surface magnetic field strength, emis- 
sion colatitude O , spectral index a , and minimum, e m i n and 
maximum, e max , energy of the primary photon spectrum. Cas- 
cade simulations have been run for the cases where one mode of 
splitting operates (the -L— >|||| mode), three modes of splitting 
operate (the _L— >||||, _L^_L_L and |— >J_|| modes permitted by 
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QED), and where splitting is turned off completely ("no split- 
ting" ) . Figure || shows the cascade pair yield as a function of 
magnetic field strength and maximum primary photon energy; 
since the total number of photons and the pair yield depend on 
the value of £ m in , we hold £ m in = 1 constant for all calcula- 
tions. 
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Fig. 6. — The cascade pair yield (number of pairs per injected pho- 
ton) as a function of surface magnetic field strength, Bo , in units of the 
critical field, B CI , and maximum primary photon energy, E max , in units 
of m e c 2 . Different line types refer to the cases where no splitting, only 
one mode ( _L— +|||| ) or three modes of splitting are allowed. The curves for 
the "no splitting" case behave similarly for e m ax — 10 4 , £max = 10 6 and 
£max = 10 s cases, so only the curves for the £ max = 10 2 and e max = 10 6 
cases are shown for clarity. Switching on only one mode of splitting leads 
to only minimal reductions in the pair yield for e max > 10 4 , as explained 
in the text. Photons are assumed to be emitted from the stellar surface. 

The curves for the "no splitting" case, shown only for e max = 
10 2 and e ma x = 10 6 to avoid confusion, are of course not physi- 
cal but are included to illustrate the effect of one splitting mode 
on the pair yield. When e max is below the pair escape energy 
(cf. Section 2.1.3 ), as it is for the £ max = 10 2 and £ max = 10 4 
curves at lower field strengths, the pair yield is considerably 
lower than when £ max is well above pair escape energy (as is 
the case for the £ max = 10 6 and £ max = 10 8 curves), for which 
pair yields at the lower field strengths can be quite large due 
to the pair multiplicatio n ef fect through successive generations, 
as described in Section 2.4. Above a field strength of 0.1i? cr , 
in £max 10 cases the pair yields drop because the pairs are 
produced in low Landau states at or near threshold, suppressing 
synchrotron emission and thus the pair multiplication effect. In 
the absence of photon splitting, the pair yield increases slowly 
with increasing Bq > 0.1i? cl , because the number of generations 
of the cascade can increase. The rapid rise of the £ max = 10 2 
case when Bq <^ 0.2B cr is due primarily to the escape energy for 

pair production £c S T e+e being close to £ max so that the opti- 
cal depth is of the order of unity or less; the pair yield must then 
necessarily be a strong function of Bq by virtue of the expo- 
nential asymptotic form (see Daugherty & Harding 1983) for the 
pair creation rate. Since the pair escape energy saturates when 
Bq <; 0.5-B cr at approximately the threshold value of 2/ sin 9 



(see Figure |l|), the pair yield should asymptote to a constant 
value at high field strengths. This constant is just the frac- 
tion of primary photons above pair threshold: when £ max 3> 1 , 
this fraction is approximately (sin9/2) Q ~ 1 , which evaluates to 
~ 0.15 for a polar cap size of = 5° and a = 1.6 . Reducing 
£max clearly diminishes the proportion of photons above pair 
threshold. 

In the presence of photon splitting, above Bq ~ 0.5£> cr , the 
cascade is limited to only one pair generation in each polariza- 
tion mode, amounting to one pair per || mode photon injected 
above the pair escape energy plus two pairs per _L mode photon 
injected above the splitting escape energy (see Figure [|). The re- 
sults clearly indicate that photon splitting in only one mode does 
not significantly suppress the pair yield. As long as one photon 
polarization mode does not undergo splitting, there is always 
a channel for pair production. In fact, photon splitting in one 
mode can even increase the pair yield, as is evident in Figure [}] 
for field strengths between 0.2£? cr and 0.7_B cr , for £ max > 10 4 . 
Consequently, the putative radio quiescence boundaries in Fig. ^| 
are not borders to domains of significant pair suppression if pho- 
tons of || polarization cannot split. While splittings -L— >|||| can 
compete effectively with pair creation, the two produced || pho- 
tons can only create pairs if splitting of || photons is forbidden 
by kinematic selection rules. In such a case, pair creation is 
prolific unless the second generation || photons are around or 

below the escape energy for pair production, £cs^ e e , 



the 



II y c 

original _L photon is of energy <; 2£e SC 



Hence, in the en- 



ergy range £o SC 



£ £ £ 2eLT e , -L-HIH splittings truly 
do suppress pair creation; otherwise they merely postpone it 
a generation. By comparing the number of primary photons 

II > e + e - 

in this range with that above £c S c , it is quickly estimated 
that the pair yield is enhanced by a factor of the order of 2 2 ~ Q 
when splitting is switched on, i.e. generally a minimal change. 
This is precisely the behaviour observed in Figure ^ for all but 
the £ max = 10 2 case. Effectively, the doubling of the number of 
photons by splittings compensates almost exactly for the sup- 



pression of pair creation below 2e, 



. Therefore, when only 



_L mode photons can split, pair creation is inhibited only when 

r> || — ►e + e — 
£max ^ ^£csc 

The situation changes dramatically if three modes of photon 
splitting (those allowed by QED) are operating. As shown in 
Figure g for the case of £ max = 10 2 and £ max = 10 4 (the other 
cases not shown behave similarly to this case), the pair yield 
drops dramatically above B ~ 0.4i? cr (depending on 9 ) when 
photons of both polarization modes are allowed to split. Pure 
photon splitting cascades can now take place, and typically ex- 
ceed ten photon generations. In fields above B CI very few pairs 
are created because the photon splitting cascade degrades the 
photon energies below the pair escape energy before the cascade 
reaches a height where pair production can take over. Figure ^ 
shows the effect of varying the primary photon injection co- 
latitude on the cascade pair yield, as a function of magnetic 
field strength for cascades where three modes of splitting are 
allowed. The field strength at which the pair yield begins to 
drop increases with decreasing colatitude, due to the increas- 
ing field line radius of curvature. This effect was seen in the 
radio quie scen ce boundary, where Bq cx 9 -4 / 15 was computed 
in Section iA by equating pair and splitting escape energies. A 
similar boundary may be predicted by equating the polarization- 
averaged escape energies for pair creation and photon splitting, 
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with the result being a virtually identical dependence of Bq 
on , but with a constant of proportionality that is a factor 
of about 1.7 lower than that in Eq. [0. Such a polarization- 
averaged boundary of radio quiescence reproduces the ( Bq , Q ) 
phase space corresponding to the precipitous declines observed 
in Figure ^. Hence it can be concluded that the simplified cri- 
terion for pair suppression by splitting is fairly accurate if all 
three QED -permitted modes of splitting are operating. 
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Fig. 7. — The cascade pair yield (number of pairs per injected photon) 
as a function of surface magnetic field strength, Bo , in units of the crit- 
ical field, B CI , for different primary photon injection colatitudes, , for 
the case where three photon splitting modes operate. Pair suppression be- 
comes substantial in near-critical and supercritical fields when pure photon 
splitting cascades operate. Fluctuation in the curves at the highest field 
strengths is numerical noise due to low counts. Again, surface emission of 
photons is assumed. 

As a more complete representation of the relevant phase 
space, contour plots of the pair yield as a function of magnetic 
field and pulsar period are shown in Figure |^ for the cases of 
one and three splitting modes. The period is derived from the 
surface injection colatitude assuming emission at the rim of a 
polar cap of standard size as given in Eq. (^). In the left 
panel of Figure |[ which shows the case for one splitting mode, 
the most dramatic feature is the sharp fall-off of pair yield at 
large pulsar periods. This effect is purely a consequence of the 
polar magnetic field geometry: at large periods (small injection 
colatitudes) the radius of curvature of the field lines is larger 
and the pair production attenuation length can exceed a stellar 
radius, suppressing pair cascades. This defines the "death line" 
for radio pulsars. There is a weak dependence of the "death 
line" on magnetic field strength for B J> 0.1 B cr because the 
pair attenuation length saturates due to threshold pair produc- 
tion. With only one splitting mode in operation, the pair yield 
has a relatively weak variation with magnetic field. However, 
there is a maximum in pair yield at fields B <J 0.1 -B cr and short 
periods, because pairs are produced in high Landau states. 

In the right panel of Figure which shows the case where 
three splitting modes are operating, the most dramatic feature is 
a sharp fall-off of pair yield at magnetic fields around B ~ B cl . 



This feature is the radi o qu iescence line due to photon splitting 
as described in Section 3J The period dependence of this line 
is consistent with that inferred from the escape energy calcula- 
tion, i.e. B oc e~ 4 / 15 oc p- 2 / 15 . The radio "death line" is 
also apparent as a decline in pair yield at large periods. 

4. DISCUSSION 
4.1. Radio quiescence? 

The results of the previous section show that if three modes 
of splitting operate at high field strengths, then there is nearly 
complete pair suppression at high field strengths and the approx- 
imate condition for radio quiescence proposed in Section 3.1 is 
valid. If only one mode of splitting (i.e. -L— >|||| ) operates, then 
some other effect must act to reduce the pair yield, otherwise, 
pair creation in cascades remains prolific at near-critical and 
supercritical fields, and there is no reason to expect a regime 
of radio quiescence in the upper portion of the P - P diagram. 
An answer to the question of how many modes of photon split- 
ting operate in highly dispersive magnetic fields is not presently 
available and will require detailed investigation into the behav- 
ior of the photon splitting rate in this regime. 

We have discussed several ways in which a suppression of pair 
creation at high B can be present if only one splitting mode op- 
erates. First, if the effective maximum energy e max of primary 
photons satisfies e max <J 2e^^ e e , then the pair yield must 
be dramatically reduced whenever such photons are mostly of 
_!_ polarization. However, there is no reason a priori to believe 
that such low values of e max should arise when B J> B cl - . For 
curvature radiation, the spectrum is only dependent on Bq via 
acceleration properties. For resonant Compton scattering, the 
spectral dependence on Bq is more complicated and needs to be 
the subject of future investigation. While the maximum electron 
energy declines somewhat with Bq (e.g. Harding & Muslimov 
1998; see the discussion just below), the primary photon spec- 
trum is flat up to ~ B 2 /e s (for soft target photon energy e s and 
drops steeply thereafter: Dermer 1990; Baring 1994). Hence, 
properties of the primary spectrum are unlikely to uniformly 
produce suppression of pair creation at high field strengths. 

A second possibility that is effective at diminishing, but not 
completely suppressing, the number of pairs when e max 3> 
£ ic~* e e i s ground state pair creation. It inhibits successive 
generations of pair cascading and becomes significant in local 
fields B ;> 2 x 10 12 Gauss. This value is below the spin-down 
fields of the most magnetic radio pulsars, so that evidently 
ground state pair creation cannot be primarily responsible for 
imposing radio quiescence. Yet this physical effect could be a 
contributory factor, noting that the extant pulsar population 
is somewhat thinner in the P-P diagram than might be ex- 
pected at Bq ;> 10 13 Gauss. A related possibility (as discussed 
in Section 2.3) is that positronium formation may inhibit radio 
emission without reductions in pair creation. But even if sta- 
ble positronium formation does occur in pulsar magnetospheres, 
it cannot account for observed radio quiescence in pulsars be- 
cause the predicted onset of positronium formation occurs at 
field strengths ~ 0.1£? cr , in the middle of the radio pulsar pop- 
ulation. Indeed, Usov & Melrose (1996) invoke bound-state pair 
creation in their model for gamma-ray emitting radio pulsars. 

Finally, we have argued that suppression of pair creation by 
the conversion of || polarization photons that are the product of 
-L— HIH splittings back into the _L state by the resonant Comp- 
ton scattering process, before they can pair produce, will not 
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Fig. 8. — Contour plots of cascade pair yield (number of pairs per injected photon) as a function of surface magnetic field strength, Bo , in units of the 
critical field, B CI , and pulsar period for the cases where (left panel) only one mode of photon splitting is allowed, and (right panel) where three photon 
splitting modes are permitted. Contour levels are equally spaced logarithmic (base ten) intervals. 



be effective in normal pulsar magnetospheres where the particle 
densities are low. 

While these theoretical issues leave the question of radio qui- 
escence quite open, several observational issues are pertinent to 
this discussion. These relate to whether or not one expects to 
observe radio pulsars with high P , which divides into issues 
of intrinsic population densities in the P - P diagram and ob- 
servational selection effects. There are several obvious natural 
biases influencing pulsar observability. First, due to the expec- 
tation that pulsars with higher spin-down power (i.e. Bq/P 4 ) 
are probably more luminous, a property observed in both the 
X-ray (Becker & Triimper, 1997) and gamma-ray (e.g. Zhang 
& Harding 2000) bands, one might anticipate that high-field ra- 
dio pulsars are more luminous than their less-magnetized coun- 
terparts. There is evidence that this trend is present in the 
observed population: when luminosity distributions for Prince- 
ton Pulsar Catalog members are binned in ranges of Bq , radio 
luminosities L ra d increase with spin-down field up to around 
2 x 10 13 Gauss, beyond which the trend appears to reverse and 
£rad perhaps begins to decline with B . While this reversal 
is suggestive of an onset of radio-quiescence at high fields, it is 
presently statistically insignificant. The addition of the Parkes 
Multi-Beam population (for which distances have not yet been 
established from the dispersion measures) will improve statis- 
tics, but perhaps still leave the evidence for a luminosity decline 
inconclusive. Uncertainty in source distances pervades estimates 
of pulsar luminosity, and since there is a correlation between 
distance out of the galactic plane and pulsar age (and therefore 
period), we caution against over-interpretation of L ra d trends 
for high Bq ■ 

The next factor is the age distribution in the P-P diagram. 
Clearly, evolution (with or without field decay) guarantees a 
denser population at long periods, and the speed at which the 
period evolves is a rapidly increasing function of Bo ■ The popu- 
lation density would be constant along diagonal lines of constant 
age P/(2P) for a uniform birth rates vs. field distribution. Un- 



less there are many fewer pulsars born with high field strengths, 
the smaller average age of high field pulsars to the left of the 
conventional death lines cannot account for the dearth of radio 
pulsars at these fields. Concomitantly, the expected clustering 
of pulsars near the death line at subcritical fields is not realized; 
the beaming of radio emission is obviously a contributor to this 
effect. The reduction of the polar cap size with period should 
reduce the solid angle ( oc P^ 1 ) of the cone of emission to more 
or less compensate the age clustering effect. Such a property 
was invoked by Young, Manchester & Johnston (1999) to ar- 
gue that the 8.5 second Parkes pulsar PSR J2144-3933 was not 
isolated in its existence, but rather representative of a large, un- 
seen population that challenges conventional theories explaining 
the death line. Yet the observed thinning of the population near 
this boundary argues that another factor is influential in deter- 
mining the period distribution; a reduction in luminosity due to 
the onset of pair quenching may provide a partial explanation 
of this dilution. Notwithstanding, the beaming phenomenon 
impacts the P distribution and not P phase space. 

There are two main observational selection effects, the first 
being an intrinsically greater sensitivity at longer periods and 
lower dispersion measures, basically due to interstellar disper- 
sion and scintillation effects, radiometer noise, and pulse shape 
and Fourier analysis properties (e.g. see Cordes and Chernoff 
1997 for a discussion of sensitivities). The second true selec- 
tion effect, less evident in the literature, is that pulsar surveys 
tend to filter out baseline fluctuations on long (i.e. J> 5 second) 
timescales, thereby selecting against supersecond period pulsars 
(Cordes, private communication). This bias is partly driven by 
past needs to focus on the appropriate range of periods of con- 
ventional pulsars, and will be counterbalanced by the inc reasing 
scientific interest in magnetars, discussed in Section Els below. 
In summation, in the light of all these observational consider- 
ations, without performing an extensive population statistics 
analysis, it is fairly safe to argue that radio pulsars with P 
in excess of 10 -11 are either rare or non-existent. The confir- 
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mation of a possible radi o pulsar counterpart to SGR 1900+14 
(discussed in Section 4.3), or otherwise, shall play a prominent 
role in resolving this issue. 

4.2. Particle Acceleration Locales 

Since the strength of the local magnetic field is one of the 
key influences on the location of the high-field death line, the 
site of the high energy emission in the pulsar magnetospherc is 
critically important. The initial angle at which the high-energy 
photons are emitted relative to the local magnetic field has also 
been shown in the previous subsection to be of critical impor- 
tance (see Figure |^) . Both of these factors are determined by the 
nature of the particle acceleration above the polar cap. Several 
types of models have studied pulsar acceleration due to charge 
deficits at different locations in the magnetosphere. Polar cap 
models consider the formation of a parallel electric field in the 
open field region near the magnetic poles, while outer gap mod- 
els consider acceleration in the outer magnetosphere, near the 
null charge surface (see Mestel 1998 for the most recent and 
comprehensive review of pulsar electrodynamics). 

Polar cap models for pulsar high-energy emission are all based 
on the idea, dating from the earliest pulsar models of Sturrock 
(1971) and Ruderman & Sutherland (1975; hereafter RS75), of 
particle acceleration and radiation near the neutron star surface 
at the magnetic poles. Within this broad class, there is a large 
variation, with the primary division being whether or not there 
is free emission of particles from the neutron star surface. This 
question hinges on whether the surface temperature T of the 
neutron star (many of which have now been measured in the 
range T ~ 10 5 — 10 6 K; Becker & Triimper 1997) exceeds the 
ion, T; and electron, T c , thermal emission temperatures. If 
T < TJ , a vacuum gap will develop at the surface, due to the 
trapping of ions in the neutron star crust (RS75, Usov & Melrose 
1995). In this case, the particle acceleration and radiation will 
take place very near the neutron star surface. If T > T e , free 
emission of particles of either sign of charge will occur. The flow 
of particles is then limited only by space charge (SCLF case), 
and an accelerating potential will develop (Arons & Scharle- 
mann 1979; Muslimov & Tsygan 1992) due to an inability of 
the particle flow all along each open field line to supply the 
corotation charge (required to short out Eu ). In space charge- 
limited flow models, the accelerating Eu is screened at a height 
where the particles radiate 7 -rays that produce pairs. This so- 
called pair formation front (PFF) (e.g. Arons 1983, Harding & 
Muslimov 1998) can occur at high altitudes above the polar cap. 

Zhang & Harding (2000b) propose that, if both photon po- 
larization modes can undergo splitting at high fields, the radio 
quiescence line will be much lower for pulsars with vacuum gaps 
("anti-pulsars" in the language of RS75), than for pulsars with 
SCLF gaps ("pulsars"). In the case of vacuum gaps the high 
energy radiation occurs near the neutron star surface, where 
the local fields are high, whereas in the case of SCLF gaps, the 
high energy emission (and subsequent pair creation) may occur 
at high altitudes, especially when photon splitting prevents pair 
creation and thus a PFF near the surface. The radio quiescence 
line for "anti-pulsars" will be the photon splitting death line for 
surface emission, while the radio quiescence line for "pulsars" 
will be located around Bq = 2 x 10 14 G (for a surface temper- 
ature of lO 6 -^), above which thermionic emission of electrons 
from the surface is no longer possible and "pulsars" must have 
vacuum gaps. This could account for the observed existence of 
radio-loud pulsars such as PSR J1814 and radio-quiet pulsars 
such as AXP IE 2259 at the same field strength and in close 



proximity in the P - P diagram. 

4.3. Magnetars 

The importance of a boundary for radio quiescence to the 
study of the radio pulsar population is obvious. Yet the absence 
of radio pulsars in the high- Bq region of phase space has become 
of even greater importance with the mounting observational ev- 
idence at X-ray and 7-ray energies for the existence of neutron 
stars with such large fields, which a priori are not discriminated 
against (for a given age ) by known radio selection effects (as 
discussed in Section 4.1). This evidence includes the detection 
of spin-down in the growing number of anomalous X-ray pul- 
sars (AXPs). These sources have been known to exist for over 
a decade (e.g. see the biographical summary of Mereghetti & 
Stella 1995), though the identification of them as being anoma- 
lous was forged slowly with the accumulation of observational 
data. The designation "anomalous," coined by van Paradijs, 
Taam & van den Heuvel (1995), was founded in their periods 
(6-12 seconds) being relatively short compared with typical ac- 
creting X-ray binary systems, combined with their unusually 
steep X-ray spectra and monotonic increases in periods: short 
period binary X-ray pulsars usually exhibit P < , i.e. spin-up. 
A list of AXPs with measured P and P is given in Table 1, 
an adaptation of that in Gotthelf & Vasisht (1998). Further- 
more, the association of two of these pulsars (IE 1841-045 and 
IE 2259+586) with supernova remnants (which are typically 
much younger than X-ray binaries) has shifted the focus from 
accretion torques, as championed by Mereghetti & Stella (1995) 
and van Paradijs, Taam & van den Heuvel (1995), to electro- 
magnetic dipole torques (e.g. Vasisht & Gotthelf 1997) as the 
origin of the spin-down. In fact, Vasisht & Gotthelf (1997) argue 
that it is difficult for accretion torques to spin a pulsar down to 
P ~ 10 seconds on a 10 3 year timescale unless the pulsar was 
born a slow rotator. 

The electromagnetic dipole interpretation for the spin-down 
immediately implies immense supercritical ( Bq > 4.41 x 
10 13 Gauss) fields in these sources: inferred values are given in 
Table 1. The motivation for such a perspective has been dramat- 
ically enhanced by the recent detection of similar periods and 
period derivatives of soft gamma repeaters (SGRs), the neutron 
star sources with transient outbursts of soft gamma-ray emis- 
sion possessing quasi-thermal spectra. Duncan & Thompson 
(1992) postulated that neutron stars with supercritical magne- 
tizations, magnetars, were responsible for SGR activity based 
on the observed 8 second periodicity (Mazets et al. 1981) of the 
5th March 1979 outburst from SGR 0525-66, combined with it's 
strong directional association (Helfand & Long 1979; Cline et al. 
1982) with the supernova remnant N49 in the Large Magellanic 
Cloud. No spin-down P was measurable for this source since 
the data that unambiguously exhibited periodicity spanned a 
single range of around two minutes. Hence the surface field of 
Bq ~ 6 x 10 14 Gauss inferred by Duncan & Thompson (1992) 
for SGR 0525-66 was purely circumstantial. Nevertheless, their 
hypothesis was dramatically bolstered by the detection of a su- 
persecond periodicity coupled with a high rate of spin-down in 
the quiescent counterparts of SGR 1806-20 (Kouveliotou, et al. 
1998) and of SGR1900+14 (Kouveliotou et al. et al. 1999); 
parameters for these sources are also listed in Table 1. A giant 
burst, very similar to that observed from SGR 0525-66, was seen 
from SGR1900+14, exhibiting the same 5 s period as was de- 
tected in the quiescent emission (Hurley et al. 1999), strength- 
ening the magnetar identification for SGRs. While the similar- 
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TABLE 1 

Spin-Down Parameters 61 for Anomalous X-ray Pulsars and Soft Gamma Repeaters 









P 




P 




r 


P/2P 




B b 


Pulsar 


SNR 


Ref. 


(sec) 


(sec 


sec - 






(G 


auss) 


SGR 1806-20 


G10.0-0.3 


c 


7.47 


8.3 


x 10 


-li 


I 


x 10 3 


1.6 


x 10 15 


SGR 1900+14 




d 


5.16 


1.1 


x 10 


-10 


7.8 


x 10 2 


1.5 


x 10 15 


IE 1841-045 


Kes 73 


e 


11.77 


4.7 


x 10 


-11 


4.0 


x 10 3 


1.5 


x 10 15 


IE 1048-5937 




f 


6.45 


2.2 


x 10 


-1 1 


4.6 


x 10 3 


7.6 


x 10 14 


4U 0142+615 




g 


8.69 


2.3 


x 10 


-12 


6.0 


x 10 4 


2.9 


x 10 14 


IE 2259+586 


CTB 109 


g 


6.98 


7.3 


x 10 


-13 


1.5 


x 10 5 


1.4 


x 10 14 


1RXS J170849. 0-400910 




h 


11 


2.3 


x 10 


-1 1 


8.0 


x 10 3 


1.1 


x 10 15 



Note. — (a) Only sources with measured P are included in the Table, (b) The magnetic field esti- 
mates are obtained using the choice for the dipole spin-down energy loss rate B = 6.4 x 10 19 (P P) 1 / 2 
Gauss adopted by Usov & Melrose (1995), derived in, for example, Shapiro & Teukolsky (1983). Ref- 
erences for observed periods and period derivatives are (c) Kouveliotou et al. (1998), (d) Hurley et al. 
(1999) for the period and Kouveliotou et al. (1999) for P for SGR 1900+14, (e) Vasisht and Gotthelf 
(1997), (f) Oosterbroek et al. (1998; here the estimate obtained from their Table 3 using the maximum 
time baseline is used), (g) Mereghetti & Stella (1995), (h) Israel et al. (1999). 



ity of P and P and the existence of remnant identifications for 
both AXPs and SGRs may provide moderately compelling moti- 
vation for classifying the two types of sources together, given the 
vastly different emission properties between AXPs and SGRs, 
their commonality may be confined just to magnetars being the 
sites for their activity. 

A property mutual to anomalous X-ray pulsars and soft 
gamma repeaters is that they are radio quiet pulsars (if one 
excepts the unconfirmed detection of SGR 1900+14 by Shitov 
1999), an intriguing fact given their high magnetic fields and 
the focus of this paper. The spin-down parameters for these 
sources (see Table 1), are used to derive their positions on the 
P-P diagram in Figure g. The AXPs and SGRs all lie in 
the extreme upper right of the P-P phase space, above our 
derived #kB radio quiescence line. However, one of the radio pul- 
sars recently discovered in the Parkes Multi-beam Survey, PSR 
J1814-1744, lies very close to the AXP CTB 109, such that 
a single radio quiescence line cannot separate the radio pulsar 
and AXP populations. This proximity coupled with the fact 
that PSR J1814-1744 has not been detected in X-rays (Pivo- 
varoff, Kaspi & Camilo 2000) strongly suggests that a quantity 
other than P and B a has a profound influence on the prop- 
erties of highly-magnetized radio pulsars and AXPs. Note that 
there are several effects which could cause a fuzziness in both 
the radio quiescence boundary and the positions of t he pulsars 
in the P-P diagram. As mentioned in Section 3.1, the effec- 
tive quiescence boundary for a particular source depends on the 
altitude of emission. The quiescence boundary applicable for 
high-B radio pulsars, whose high-energy em issio n may occur at 
some altitude above the surface (see Section 12), may be higher 
than the quiescence boundary applicable for AXPs, whose emis- 
sion is thought to be predominantly of thermal origin from the 
hot neutron star surface. Another possibility, pertinent if the 
radio detection of SGR 1900+14 is confirmed, is that SGRs 
are "pulsars" with active accelerators allowing pair production 
and radio emission at high altitudes, as discussed by Zhang & 



Harding (2000b, see Section 4.2), while AXPs are radio-quiet 
"anti-pulsars" . 

Furthermore, the actual magnetic fields of the AXP and SGR 
sources are not likely to be those derived from pure dipole spin- 
down. If rapid field decay is powering the luminous emission 
from these sources (Thompson & Duncan 1996), then it is likely 
that transient higher multipoles exist near the surface. Since 
higher multipole radiation contributes relatively less to the spin- 
down torque, multipole fields significantly stronger than those 
derived assuming a pure dipole could be present at the sur- 
face without affecting the spin-down. This would effectively 
lower the radio quiescence boundary for the AXPs, both due to 
the increase in field strength and decrease in radius of curva- 
ture. Measurement of braking indices larger than 3 (the dipole 
value) in these sources may reveal the presence of higher mul- 
tipoles. In fact, a marginally significant detection of v in the 
AXP 1RXS J170849.0-400910 (Kaspi, Chakrabarty & Stein- 
berger 1999) would give a braking index in excess of 100 ! Field 
decay alone would not change estimates of the surface field, but 
would decrease the age of the pulsar relative to the dipole value 
(Colpi, Geppart & Page 1999). 

On the other hand, the spin-down of some magnetars could 
be influenced by powerful particle winds, in this case lowering 
the estimated surface field strength and increasing the charac- 
teristic age (Harding, Contopoulos & Kazanas 1999, Thompson 
et al. 1999) because the wind increases the spin-down rate by 
distorting the dipole field inside the light cylinder. This is more 
likely to be a factor in SGRs, where there is evidence for par- 
ticle emission associated with bursts (Frail, Vasisht & Kulkarni 
1997; Frail, Kulkarni & Bloom 1999). Wind luminosities in ex- 
cess of ~ 10 36 ergs _1 in SGR1806-20 would lower the derived 
surface field strength below 10 14 G. An episodic wind, which 
is more likely, would decrease the estimated surface field to a 
lesser degree. Possibly for this reason, the dipole surface fields 
of the SGRs are higher than those of most of the AXPs and 
more than an order of magnitude higher than the most strongly 
magnetized radio pulsars. 
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Fig. 9. — The upper region of the P-P diagram, with filled circles again denoting the locations of members of the latest edition of the Princeton 
Pulsar Catalogue, and the open circles marking pulsars in the recent Parkes Multi-Beam survey. The Crab, Vela and PSR 1509-58 gamma-ray pulsars 
are highlighted as indicated in the lower left inset, together with the positions of five radio-quiet anomalous X-ray pulsars (open square stars) and SGRs 
1806-20 and 1900+14 (filled squares) in the upper right of the diagram; their measured P and inferred fields So are listed in Table 1. The conven- 
tional death line (as in Fig. |j) lies at periods longer than all seven of the putative magnetars. The 0kB,o = (midweight solid diagonal curve) and 
$kB,o = 10~ 3 (midweight dashed curve) depictions of the radio quiescence boundary from Fig. ^ are exhibited. The heavy-weight pairs of curves labelled 
£max = lOMeV and e max = 10 2 MeV are contours for the escape energy of photon splitting, _L— *|||| , with the solid and dashed line styles denoting 
$kB = an d #kB,o = 10 -3 , respectively. They are shown only above the potential radio quiescence boundaries, where splitting dominates pair creation 
for the X polarization, marking contours for the approximate maximum observable energy permitted in highly-magnetized pulsars; to the right of these 
the magnetosphere is transparent to 10 MeV and 100 MeV photons, respectively. 



4.4. Guides for Soft and Hard Gamma-Ray Pulsation Searches 

To close this section of discussion, it is appropriate to mention 
the implications of our analysis for gamma-ray pulsar searches, 
which are, in a sense, independent of the radio quiescence is- 
sue. If the electrons are accelerated to high Lorentz factors in 
profusion, then primary photons will extend into the EGRET 
band and beyond. Therefore, short period pulsars with high 
spin-down power ( oc Bq /P 4 ) like the Crab pulsar should be 
luminous and therefore easily detectable by future experiments 
such as GLAST. At longer periods near the magnetar regime, 
the collecte d ob servational data and the age distribution bias 
(see Section 4.1) indicate that spin-powered pulsars will be gen- 
erally faint unless they are nearby. Yet the observed hard X-ray 
luminosities of AXPs are well in excess of their spin-down lu- 
minosities so that non-rotational energy sources are indicated. 
Hence luminosity biases against long period pulsars may not be 
as pronounced for those with alternative power sources such as 
magnetic field decay, at least in the magnetar regime. The work 
of this paper clearly underlines the fact that photon attenuation 
considerations have as profound an influence on gamma-ray ob- 
servability as do issues of power in the primary radiation. 



At near-critical and supercritical fields, the attenuation by 
pair creation and photon splitting in the EGRET band is dra- 
matic, as is evident from the contours of maximum energy de- 
picted in Figure || Basically, photon transparency at a given 
energy is achieved to the right of and below (i.e. at longer peri- 
ods and lower P ) a given contour so that the positioning of the 
e csc = 100 MeV splitting escape energy contours indicates that 
it will not be easy for an experiment like GLAST (depending on 
its final design) to detect many or most of the known magnetars, 
and the new Parkes Multi-Beam sources PSR J1119-6127 and 
PSR J1726-3530. Should magnetars emit in the gamma-rays, 
they can only do so generally below the EGRET band. This 
conclusion holds even if the magnetosphcric field structure is 
non-dipolar, since escape energies are pushed to lower values by 
greater field curvature. A possible evasion of this constraint is 
that the inferred fields are significant overestimates of the true 
surface fields, as m ight occur for wind-aided spin-down in mag- 
netars (see Section 4.3), so that the escape energy contours may 



possess pathological distortions in the P-P diagram. More- 
over, such a two-dimensional phase space may be insufficient to 
account for the magnetospheric and photon attenuation prop- 
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crties. Hence high energy astrophysicists interested in pulsar 
searches at high P should extend their period range to include 
supersecond periods to maximize their potential harvest. While 
this range has historically been given low priority (e.g. see the 
search in Mattox et al. 1996), the excitement generated by AXP 
and SGRs in recent years is changing emphases in such search 
programs. 

These attenuation properties emphasize that spectroscopy in 
a variety of gamma-ray pulsars will probably provide the ability 
to discriminate between the applicability of polar cap or outer 
gap models. In the polar cap picture, the analysis of this paper 
shows that spectral cutoff energies possess a strong inverse cor- 
relation with surface field strength Bo and the polar cap size . 
For Crab-like and Vela-like pulsars such cutoffs are coupled to 
magnetic pair creation, and appear in the EGRET band. How- 
ever, for highly-magnetized pulsars such as PSR 1509-58, such 
properties are ideally explored with a medium energy gamma- 
ray experiment, and probe the action of photon splitting. Ex- 
pectations for trends of gamma-ray cutoffs in outer gap models 
are not as well studied. Yet, it appears that their trends with 
B and pulsar period should differ significantly from those for 
polar cap models, due to the inherently different physics in- 
volved. Cutoffs in the outer gap scenario represent maximum 
energies of acceleration rather than the effects of photon absorp- 
tion, and so may actually be independent of or increase with 
surface field strength, and decline with pulse period, contrary 
to the indications of the polar gap model. Such distinctive pre- 
dictions can be probed by significant population datasets such 
as those to be afforded by the GLAST mission. A related obser- 
vational diagnostic is spawned by the contention (Chen & Rud- 
erman 1993) that gamma-ray emission is not expected in outer 
gap models when the pulsar period exceeds a sizeable fraction 
of a second. Hence detection of gamma-ray pulsations from a 
source with a supersecond period would clearly favor the polar 
cap model. 

Our studies also indicate that the gamma-ray spectra of 
pulsars should produce distinctive polarization signatures for 
the polar cap model. This should enable discrimination from 
outer gap scenarios in a future era when gamma-ray polarime- 
try is possible. While the predictions of gamma-ray polar- 
izations might be similar for the continuum spectra in each 
of these two competing models, principally because their con- 
tinuum emission processes are similar, the pair creation and 
photon splitting mechanisms generate spectral cutoff energies 
that are polarization-dependent (as emphasized by HBG97), 
and therefore immediately distinguishable from outer gap model 
turnovers that have no intimate connection to QED processes 
in strong fields. In particular, when Bo ^ B CI , the dominance 
of splittings -L— +IIH will guarantee that the cutoff energy for || 
photons exceeds that for ones of _L polarization. In contrast, 
at lower fields strengths, photon splitting is removed from the 
picture and the escape energy for pair creation by || photons 
is slightly lower than that for _L ones, resulting in a domi- 
nance of _L photons near the maximum observable energy. Be- 
cause the strong dependence of the spectral high-energy cutoff 
in fields Bo J> B cr occurs exclusively when only one mode of 
splitting is operative, it may be possible to resolve by observa- 
tion the question of whether one or several modes of splitting 
operate in high fields. The splitting and pair production cutoffs 
are expected around 100 MeV in long period pulsars (slightly 
higher if the emission region is above the surface). If only the 
-L— >|||| mode operates, the splitting cutoff will occur at an en- 



ergy around a factor of 2-3 lower than the pair production cutoff. 
GLAST (depending on sensitivity below 100 MeV) or a future 
medium energy 7-ray detector would then see 100% || polariza- 
tion of the highest energy photons, i.e. those between ei^"" 

and £cs^ e e • On the other hand, the absence of this signature 
would be consistent with, but not necessarily imply, at least two 
active splitting modes. These distinctive signatures establish a 
strong case for performing gamma-ray polarimetry experiments, 
an objective that may not be that far in the future for medium- 
energy Compton gamma-ray telescopes. 

5. CONCLUSION 

We have presented a detailed study of the comparative at- 
tenuation of photons by magnetic pair production and photon 
splitting, of pair cascades and of conditions for pair suppres- 
sion in very highly magnetized pulsars. While ground-state pair 
creation and positronium formation can act at B <; 0.1i3 cr 
to reduce the number of free pairs, only photon splitting has 
the potential capability of dramatically inhibiting the creation 
of pairs, bound or free. Our quantitative study of cascade pair 
yields at different field strengths confirms the location of the 
approximate radio quiescence boundary in the P - P diagram 
proposed by Baring & Harding (1998), assuming that copious 
pair production is a requirement for pulsar radio emission. How- 
ever, it also shows that the existence of such a boundary requires 
the rate of photon splitting to be non-zero for both photon po- 
larization modes. This is not believed to be true in the low-field, 
weakly dispersive limit, where kinematic selection rules prohibit 
splitting of photons of || polarization. Whether photon split- 
ting can significantly suppress pairs in highly magnetized radio 
pulsars depends on presently unexplored physics. It therefore 
has become critically important to study the photon splitting 
process in the high-field, highly dispersive regime. Understand- 
ing the behavior of splitting in this regime will not only resolve 
the radio quiescence question, but is crucial to modeling accel- 
eration and radiation in high-field pulsars and magnetars. 

The circumstantial evidence, that the most highly magnetized 
pulsars (AXPs and SGRs) tend to be radio quiet, argues for 
some pair suppression at high field strengths. The radio pulsar 
and magnetar populations are no longer cleanly separated in P- 
P space, which indicates that another dimension of phase space 
is necessary to divide these two source populations. Delineation 
of these two source groups is complicated by the likelihood of 
alternate spin-down mechanisms and evolution issues for mag- 
netars. More systematic surveys of the overlap region of P - P 
space at X-ray and 7-ray energies, correlated with radio pulsar 
surveys are needed. We have given some guidelines for such 
high-energy searches based on the expected physics of particle 
acceleration and photon attenuation, in anticipation of the next 
generation of medium and high-energy gamma-ray experiments 
such as GLAST. Spectral and polarization observations may be 
crucial in identifying signatures of photon splitting and even 
in resolving the issue of which splitting modes are operating 
in ultra-magnetized sources, an attractive possibility for both 
astrophysicists and physicists. 
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